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Experiments have widely shown that a steady-state lamellar eutectic solidification front is destabi- 
lized on a scale much larger than the lamellar spacing by the rejection of a dilute ternary impurity 
and forms two-phase cells commonly referred to as "eutectic colonies". We extend the stability 
analysis of Datye and Langer [V. Datye and J. S. Langer, Phys. Rev. B 24, 4155 (1981)] for a 
binary eutectic to include the effect of a ternary impurity. We find that the expressions for the crit- 
ical onset velocity and morphological instability wavelength are analogous to those for the classic 
Mullins-Sekerka instability of a monophase planar interface, albeit with an effective surface tension 
that depends on the geometry of the lamellar interface and, nontrivially, on interlamellar diffusion. 
A qualitatively new aspect of this instability is the occurrence of oscillatory modes due to the in- 
terplay between the destabilizing effect of the ternary impurity and the dynamical feedback of the 
local change in lamellar spacing on the front motion. In a transient regime, these modes lead to 
the formation of large scale oscillatory microstructures for which there is recent experimental evi- 
dence in a transparent organic system. Moreover, it is shown that the eutectic front dynamics on 
a scale larger than the lamellar spacing can be formulated as an effective monophase interface free 
boundary problem with a modified Gibbs-Thomson condition that is coupled to a slow evolution 
equation for the lamellar spacing. This formulation provides additional physical insights into the 
nature of the instability and a simple means to calculate an approximate stability spectrum. Finally, 
we investigate the influence of the ternary impurity on a short wavelength oscillatory instability that 
is already present at off-eutectic compositions in binary eutectics. 



PACS: 81.30.-t, 81.30.Fb, 64.70.Dv 
I. INTRODUCTION 

The interfacial patterns that arise naturally during the 
solidification of eutectic alloys have attracted widespread 
interest for several decades from both fundamental and 
practical viewpoints. At a fundamental level, the main 
theoretical challenge lies in understanding the complex 
spatiotemporal dynamics of phase boundaries (solid- 
liquid and solid-solid) resulting from the competition of 
two thermodynamically stable solid phases growing si- 
multaneously into a metastable liquid phase. In particu- 
lar, one basic question is how to understand the nature of 
the morphological instability of the simplest spatially pe- 
riodic steady-state that gives rise to this rich dynamics. 
From a practical viewpoint, the composite microstruc- 
ture formed by lamellae or rods of these two solid phases 
growing simultaneously from the melt leads to interesting 
materials where the properties of two different solids can 
be advantageously combined. Moreover, the typical size 
of the microstructure pattern is about an order of magni- 
tude smaller than in dendritic alloys, leading to superior 
mechanical properties. Consequently, understanding the 
solidification processing conditions that yield a particular 
eutectic microstructure is a goal of direct technological 
relevance. 

Since the early investigations of eutectic alloys, it has 
been remarked that besides the fine lamellae or rods, 
there may exist cellular structures, termed colonies [Q-^). 
Their size is typically 10 to 100 times the lamellar spac- 
ing. In Fig. p] we show an experimental picture of colonics 




FIG. 1. Eutectic colonies in the transparent organic alloy 
CBr4-C2Cl6 (from Ref. Q), grown by a directional solidifica- 
tion experiment. The growth direction is from bottom to top. 
The envelope of this two-phase structure (i.e. the solidifica- 
tion front on a scale much larger than the lamellar spacing) 
closely resembles the typical monophase cells observed during 
the directional solidification of a dilute binary alloy. 
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obtained in a directional solidification experiment B . In 
their overall shape, the two-phase cells are remarkably 
similar to the monophase (i.e. single solid phase) solidi- 
fication cells formed in standard directional solidification 
of a dilute binary alloy. This analogy is further supported 
by the experimental finding that colonies appear only 
when a ternary impurity, rejected by both solid phases, 
is present In contrast, in binary eutectics the large 

scale solidification front stays planar for a range of com- 
positions around the eutectic point. This suggests that 
the mechanism of the instability is similar to the classical 
Mullins-Sekerka instability of a monophase solidification 
front p0| |. Indeed, its onset is relatively well described 
quantitatively by the constitutional supercooling crite- 
rion jnj, according to which the interface becomes un- 
stable when the ternary impurity concentration gradient 
in the liquid ahead of the interface exceeds a critical value 
set by the ratio G/v p of the temperature gradient and the 
pulling speed of the sample. However, the spatiotemporal 
character of the linear modes associated with this insta- 
bility has never been investigated. In particular, it has 
remained unclear how the Mullins-Sekerka analysis has 
to be modified to account for the composite structure 
of the interface. To answer this question, we present in 
this paper a linear stability analysis of a lamellar eutectic 
solidification front in the presence of a ternary impurity. 

It is useful to first briefly review the progress accom- 
plished to date on the related problem of binary eu- 
tectic stability (without a ternary impurity). The ap- 
proach that we shall adopt here builds on earlier work 
in this context. Theoretical developments have mostly 
focused on the lamellar morphology in thin-film geom- 
etry, as then the problem can be treated as quasi two- 
dimensional. For a stability analysis, one must first ob- 
tain a steady-state solution: a shape-preserving solidifi- 
cation front propagating at constant velocity. Studies of 
this problem |ll - |l5|] led to the insight that there exists a 
family of steady-state solutions that can be parametrized 
by the lamellar spacing A. The average undercooling AT 
of the solid-liquid interface with respect to the eutectic 
temperature depends on A, and the curve AT versus A 
presents a minimum for a certain spacing A m i n - The ex- 
perimentally observed spacings are usually close to A m i n 

A hypothesis attributed to Cahn by Jackson and Hunt 
|L5j is that the lamellae always grow normal to the enve- 
lope of the solidification front. This hypothesis seems to 
work well in practice when the surface energy anisotropy 
of the solid-liquid and solid-solid phase boundaries is 
small enough to prevent locking of the lamellae to pre- 
ferred growth directions. When used to analyze hcuris- 
tically the long- wavelength stability of a eutectic front, 
this hypothesis leads to the conclusion that the lamellar 
structure is unstable for lamellar spacings below A m i n . 
Under Cahn's hypothesis, the lamellar spacing in a con- 
cave part of the solid-liquid interface decreases as the 
interface advances. Consequently, if the average spac- 
ing falls below A m i n , the local undercooling will increase 



in such a way that thinner lamellae fall further behind 
the front, leading finally to lamella termination. On the 
other hand, for spacings larger than A m i n , the opposite 
occurs: finer lamellae grow faster than wider ones and 
the concavity of the eutectic front is smoothed out. 

This argument can determine only a lower bound for 
A. To assess stability for A > A m i n , a more involved anal- 
ysis is required. Several authors tried to adapt the linear 
stability analysis of Mullins and Sekerka for single-phase 
solidification [[To) to eutectic systems. The eutectic prob- 
lem, however, is considerably more difficult because the 
basic steady-state solution is already periodic in space. 
Moreover, the presence of mobile trijunction points be- 
tween three phases complicates enormously the stability 
calculation by ruling out a smooth sinusoidal perturba- 
tion. For this reason, early attempts to average over the 
properties of the two solid phases ||17| , or_to consider per- 



turbations with immobile trijunctions 18-2^], did not 
produce consistent results (see Ref. for a more de- 
tailed discussion). 

The most complete analytical stability analysis of a eu- 
tectic interface has been performed by Datye and Langer 
(DL) ^2|. Their calculation is a perturbation analysis 
of the Jackson-Hunt (JH) jl5| steady-state solution, us- 
ing as basic variables the coordinates of the trijunction 
points both parallel and perpendicular to the interface. 
They first calculate an approximate solution to the dif- 
fusion equation for a perturbed lamellar interface. The 
assumption of local equilibrium at the solid-liquid inter- 
face and the use of Cahn's hypothesis then allow one 
to obtain an eigenvalue problem for the linear growth 
modes and to extract the stability spectrum of the inter- 
face. In the limit where the wavelength of the perturba- 
tion is large compared to the lamellar spacing (referred 
to hereafter as the "long- wavelength limit"), a simplified 
calculation confirms JH's conclusion that lamellar spac- 
ings below A m in are unstable p3| . In addition, the DL 
analysis predicted the occurrence of an oscillatory insta- 
bility with a wavelength twice the lamellar spacing for 
sufficiently off-eutectic compositions (2A-0 instability). 

The existence of this short-wavelength instability was 
later confirmed by numerical simulations of eutectic front 
dynamics using a random walk algorith m [M[ and, more 
recently, a boundary integral approach |2l|25[. The lat- 
ter study pinpointed the existence of additional short- 
wavelength instabilities, one of which (tilt bifurcation) 
was previously known |24|,|26|], and made specific quan- 
titative predictions for the CBr4-C2Cl6 organic system 
that have been validated by a detailed comparison with 
experiments [^7| , p8| . As an additional result, which is 
relevant for the present analysis, the boundary integral 
study revealed that the stability predictions of the DL 
analysis are quite accurate for lamellar spacings close 
to A m in, and only become inaccurate for larger spacings 
where the JH description of the diffusion field breaks 
down. 

In summary, according to both theory and experiment, 
a planar lamellar eutectic front in a binary alloy is com- 
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pletely stable for compositions sufficiently close to the eu- 
tectic composition, and for lamellar spacings near A m i n . 
To understand the instability leading to colony forma- 
tion, one must therefore include a ternary impurity. To 
date, few theoretical studies of ternary systems have been 
available. Rinaldi, Sharp, and Flemings derived a gen- 
eralized constitutional supercooling criterion for ternary 
systems and used it to interpret their experiments |9|. 
McCartney, Hunt, and Jordan adapted the JH steady- 
state analysis to include impurities. They interpreted the 
formation of colonies as the result of a Mullins-Sekerka 
instability driven by the ternary impurity p9| ], but did 
not carry out a detailed stability analysis. 

In this paper, we extend Datye and Langer's linear sta- 
bility approach, based on a Jackson- Hunt approximation 
of the diffusion field, to include the effect of a ternary 
impurity. For the reasons mentioned above, we expect 
this approach to yield relatively accurate predictions for 
spacings close to A m j n , which is typically the dynamically 
selected range of interest in experiments. We obtain the 
full linear stability spectrum of the steady-state lamellar 
eutectic front growing in two dimensions. Our final result 
is quite complicated, but can be substantially simplified 
for a model alloy with a symmetric phase diagram, solid- 
ified at its eutectic composition. From the study of this 
special case, we can identify all important factors that 
determine the stability of the front. In particular, we 
find that the intcrlamellar eutectic diffusion field gives a 
stabilizing contribution with a functional form similar to 
the usual capillary term. Thus, this contribution leads 
effectively to a "renormalization" of the capillary length. 
Using this insight, we are able to reformulate the stabil- 
ity problem by treating the large scale dynamics of the 
eutectic front as an 'effective monophase interface', as 
suggested by Fig. [I], with a Gibbs-Thomson condition 
that is coupled to an equation of motion for the local 
lamellar spacing. A similar type of approach has been 
used previously to analyze the long-wavelength modes of 
cellular arrays during directional solidification of dilute 
binary alloys fl3(i| ]. More recent numerical work, how- 
ever, has shown that the modes that limit the range of 
stable cell spacings at low velocity are oscillatory and 
non-oscillatory instabilities with a wavelength equal to 
twice the cell spacing [[H] that have been observed in 
experiments [ p2| . In contrast, here, no short-wavelength 
instabilities are present near the eutectic composition for 
A close tp A m i n . Therefore, the effective interface ap- 
proach provides an accurate description of the interface 
dynamics in the limit of perturbation wavelengths much 
larger than the lamellar spacing. Moreover, it can be ex- 
tended to derive a simplified expression for the stability 
spectrum by incorporating phenomenologically the effect 
of surface tension. This formula is found to predict all 
important features of the instability and to yield reason- 
ably good quantitative predictions. 

Our calculation both confirms the expectations based 
on the analogy between a two-phase eutectic front in the 
presence of a ternary impurity and a monophase front, 



and at the same time yields a surprising result. Namely, 
we find that the expressions for the onset velocity and 
wavelength of the instability are analogous to those for 
a monophase front with a surface tension renormalizcd 
by the geometry of the lamellar front and intcrlamellar 
diffusion. Thus, as far as these quantities are concerned, 
the lamellar structure leads to quantitative differences, 
but no new qualitative features of the instability. The 
new ingredient, however, which could not have been ex- 
pected on the basis of the analogy with a monophase 
front, is that the instability is oscillatory. The origin of 
this difference is due to the additional degrees of freedom 
associated with the underlying lamellar structure of the 
interface. According to Cahn's hypothesis, the change in 
the local lamellar spacing is determined by the shape of 
the front. The spacing, in turn, is related to the local 
interface temperature. As a consequence of the inter- 
play between this effect and the instability driven by the 
impurities, long-wavelength perturbations may oscillate 
during growth or form traveling waves. There indeed 
seems to be recent experimental evidence for such large 
scale oscillatory behavior near the onset of colony forma- 
tion in a transparent organic model alloy |33| . 

We also investigate the influence of the ternary im- 
purity on the short wavelength (2A-0) oscillatory insta- 
bility that is already present in a binary eutectic. The 
main result is that this instability is enhanced by the 
impurity boundary layer, which leads to a reduced com- 
position range for stable lamellar growth even below con- 
stitutional supercooling. 

The structure of this paper is as follows. In Sec. II, we 
introduce the basic sharp-interface equations. We then 
summarize in Sec. Ill the Jackson-Hunt approach and 
apply it to calculate the steady-state solution in the pres- 
ence of a ternary impurity. In Sec. IV, we review the 
principles of the DL approach and calculate the addi- 
tional terms arising from the presence of an impurity. 
Section V is devoted to a detailed discussion of the sta- 
bility spectrum at the eutectic composition in a model 
phase diagram that is symmetric about this composition. 
In Sec. VI, we reformulate the stability problem in terms 
of an effective interface approach and derive an approx- 
imate expression for the stability spectrum for an arbi- 
trary phase diagram and material parameters. In Sec. 
VII, we discuss how the off-eutectic short-wavelength os- 
cillatory instability is affected by the ternary impurity. 
Finally, we summarize our main results in Sec. VIII. 

II. BASIC EQUATIONS 

We study the solidification of a binary eutectic alloy con- 
taining a small amount of ternary impurity. Let c denote 
the concentration (in molecules per unit volume) of one 
of the constituents of the binary eutectic, and c the con- 
centration of the ternary impurity. As we restrict our at- 
tention to small impurity concentrations, we shall assume 
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that these two quantities can be treated as independent 
variables. In other words, we assume that the phase di- 
agram of the binary eutectic is only slightly altered by 
the presence of the impurity. The two solid phases are 
denoted by a and (3; c a and cp are the concentrations lim- 
iting the eutectic plateau in the binary phase diagram, 
Ac = cp — c a , and ce is the eutectic composition. 

For coupled eutectic growth, the temperature of the 
lamellar front is close to the eutectic temperature, and 
the composition on the liquid side of the solid-liquid in- 
terface is close to the eutectic point p5| . This allows to 
introduce two further simplifications. First, we may ap- 
proximate the solidus and liquidus surfaces in the ternary 
phase diagram by planes around the eutectic point. We 
denote by m s and rh s (s = a, (3) the magnitude of the liq- 
uidus slopes along the c- and c-axis, respectively. Second, 
we will assume that c a and cp are independent of tem- 
perature and impurity concentration. As already argued 
by DL, this should only slightly affect the final results, 
because in the temperature range explored by the front 
the relative variations of the concentration jumps across 
the interfaces are negligible. On the other hand, our cal- 
culations are considerably simplified as we can relate the 
volume fraction r\ of the a solid to the composition of the 
melt far ahead of the interface, Coo, via the relation 



T 

T m 



(a) 






ct+ liquid \ / p + liquid 
a + fJ + liquid ^ -.. 



c a i] + cp(l - rj), 



(1) 



FIG. 2. Phase diagram of an idealized ternary eutectic al- 
loy, (a): binary eutectic phase diagram (5=0). T m is the 
melting temperature of the pure phases, Te the eutectic tem- 
perature, and Ce the eutectic composition, (b): cut through 
the ternary diagram along the eutectic valley (c = ce)- (c): 
liquidus and solidus surfaces in the space (c, c, T). The dashed 
line is the projection of the eutectic valley on the solidus sur- 
face. The liquidus surfaces have metastable extensions be- 
yond the eutectic valley (not shown), (d): Coexistence curves 
for a fixed temperature below Te in the (c, c)-plane. 



independently of the concentration of the ternary impu- 
rity and the interface temperature. For the impurity, we 
will work in a dilute alloy approximation where the im- 
purity concentrations on the solid and liquid sides of the 
interface are related in equilibrium by 



c s = k s c L , s — a, (3. 



(2) 



The resulting phase diagram in the space (c, c, T) is 
sketched in Fig. ||. The line of intersection of the two 
liquidus surfaces is usually termed the eutectic valley or 
eutectic trough. Along this line, c, c, and T are related 
by the equations 



ce 



T-T f 



mp 



mp 



Mc 



M = 



m a mp + m a rhp 
m a + mp 



(3) 



(4) 



(5) 



Equation (g|) defines a monovariant line in the ternary 
phase diagram where, except for the special case fh a = 
rhp, the eutectic composition is shifted with respect to 
the binary eutectic point in the presence of a ternary im- 
purity. M is the liquidus slope along the eutectic valley. 
A liquid satisfying Eqs. ([}]) and (||) can be in simultane- 
ous equilibrium with two solids. 



In a typical directional solidification experiment, the 
sample is pulled in a temperature gradient G with a con- 
stant pulling speed v p . We assume that heat diffusion is 
much faster than chemical diffusion and that the ther- 
mal conductivities of solid and liquid are of comparable 
magnitude. Under this set of assumptions (commonly re- 
ferred to as the frozen temperature approximation) , the 
temperature is given by 



T(z) =T E + Gz, 



(6) 



where we have chosen the origin of the z axis at the 
eutectic temperature. 

In the absence of convection, the growth of the solid 
is limited by chemical diffusion of the constituents. We 
assume zero diffusivity in the solid (one-sided model). In 
the liquid, the diffusion equations in the laboratory frame 
(moving with velocity v p with respect to the sample) are 



1 dc 
Ddt 



d,c- 



1 dc 2 
— — = -o z c 
Ddt l 



V 2 < 



V 2 c, 



(7a) 
(7b) 



with the diffusion lengths I = 2D/v p and I = 2D/v p , D 
and D being the diffusivities of the eutectic components 
and the ternary impurity, respectively. 

For lamellar eutectic growth in thin-film geometry, the 
problem is essentially two-dimensional. Let the position 
of the solid-liquid interface be described by the curve 
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£(x,t). During the phase transformation, impurities and 
the minor component of the growing solid are rejected 
into the liquid. The condition of mass conservation im- 
plies that at the interface 

— Dd n c — v n [c(x, C) — c a ] (ah — interface), (8a) 

—Dd n c = v n [c(x, C) — cp] (PL — interface), (8b) 

—Dd n c = v n (l — k a )c(x, () (aL — interface), (8c) 

—Dd n c = v n (l — kp)c(x, C) (/3L — interface), (8d) 

where v n and d n denote the normal velocity of the inter- 
face and the derivative normal to the interface, respec- 
tively. 

We are interested in a regime of relatively low solid- 
ification velocity where the solid-liquid interface can be 
considered to be in local equilibrium, in which case the 
temperature and concentration fields at the interface are 
related by the Gibbs-Thomson conditions 



T = Te - m a (cL - ce) - fh a CL 

— T a K[(] (aL-interface) , 
T = T E + mp(c L - c E ) - rhpc L 

- TpK[Q (/^interface). 

Here, K[Q is the local interface curvature, and 
T s =T e Jsl/L s , s = a,f3, 



(9a) 
(9b) 

(10) 



are the Gibbs-Thomson constants, with 7 S £ and L s de- 
noting, respectively, the liquid-solid surface tensions and 
the latent heats of the two solids at the eutectic point. 

Finally, local equilibrium also implies that at the tri- 
junction points, where the three phases are in contact, 
the angles between the three interfaces are fixed by the 
balance of surface tension forces, which for isotropic in- 
terfaces yields the two conditions 



j aL sin a + jp L sin Op = j a p , 



7 Qi COS 8 a = Jp L COS Op, 



(11a) 



(lib) 



where the definition of the angles a and Op is illustrated 
in Fig. |. 



III. STEADY-STATE SOLUTION 

For a binary eutectic, the steady-state problem has been 
treated by Jackson and Hunt Their method has 

been extended to ternary systems by McCartney, Hunt, 
and Jordan (29). We need the steady-state solution as 
the starting point for our stability analysis. We will 
only summarize here the essential steps of the calcula- 
tion; more details can be found in Refs. [Jl5|,p2 29 1. 



A typical configuration for a lamellar eutectic growing 
at constant solidification speed v p is sketched in Fig. 
Alternating lamellae of the a and f3 phase are regularly 



liquid 



a 



X 



a 



Q+vi)X x 

FIG. 3. Sketch of a steady-state array of lamellae growing 
parallel to the z-axis. The lamella pairs are numbered by the 
integer j; rj is the volume fraction of a-phase, and z is the 
z-coordinate of the trijunction points. 



spaced. The width of one lamella pair, A, is usually of 
the order of 10 /im. We must find the interface shape for 
which the diffusion equations and the Gibbs-Thomson 
conditions at the interface are both satisfied. The JH 
approximation starts by solving the diffusion problem for 
a flat interface (consisting of lamellae with a = Op = 0). 
This can be achieved by expanding the diffusion field in 
Fourier modes along the x axis: 



c(x, z) 



^ B n exp[zg„x - q n (z - z)}, (12) 



with q n = 27m/A, and z the z coordinate of the trijunc- 
tion points. The constants q n are obtained by substitut- 
ing the above sum into the diffusion equation. Since time 
derivatives are zero in steady-state, this yields at once: 



q n = l/l + VWTqJ. 



(13) 



Inserting the expansion (12) into the mass conservation 
condition allows us to determine all Fourier coefficients 
except for Bq. The average undercooling for each lamella 
can then be calculated using the Gibbs-Thomson condi- 
tions. Finally, the condition that the two phases must 
grow at equal undercooling determines Bq and the av- 
erage undercooling AT of the solid-liquid interface as a 
function of A and 77. 

For a eutectic with a ternary impurity, we must in 
addition treat the diffusion of the ternary impurity in the 
liquid phase. For this purpose, we use the same Fourier 
expansion as above 



c(x, z) 



E 

71= — OO 



B n exp[iq n x - q n (z - z)], (14) 



where Coo is the impurity concentration far from the in- 
terface, and the constants q n are equivalent to q n with 
I replaced by I. This expansion is inserted into the con- 
dition for impurity conservation at the interface. To ex- 
tract an equation for the Fourier coefficients B n , both 
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sides of the equation are then multiplied by exp(— iq m x) 
and integrated over x from to A. The result for the 
coefficient Bo, which gives the magnitude of the overall 
diffusion boundary layer, can be written in the form 



{ T7 1 



with an effective partition coefficient 
k E = r)k a + (1 - rfjkp. 
For n 0, we obtain 

q n B n = | - ^±^j B n + l{kp- k a ) 

J- e -iqnVX_ 1 1 

X < Coo h (57 - tt)S„ 



(15) 



(16) 



-iXq n 

E — -»> 



iX(q n - q m ) 



(17) 



We are interested in a growth regime where the diffu- 
sion length is much larger that the lamellar spacing, and 
the Peclet number, Pe = X/l -C 1. Consequently, for 
0, we have q n ^> l/l, and all terms containing the 
coefficients B n on the RHS of Eq. @ can be neglected, 
which yields 



4S 



B n 



(hp - *«,) e-^" A /2 sin^A/2) 
Xq n q n l 



(18) 



Note that, if we want to go beyond this approximation, 
the problem becomes considerably more involved, be- 
cause then all the Fourier coefficients are coupled. 

We proceed now as JH by calculating the average im- 
purity concentration in front of each phase, that is, by 
evaluating 



1 f r,x 

(c)a = — T / c(x,z)d,X, 

VX Jo 



1 



A 



(1 - r))X J vX 
at z = z. The results are 



c(x, z)dx 



c 2A 
kE lr) 



2A 



k E 1(1 -v) 



P(i])coo{k[3 - k a 



with 



oo _^ 

p (v) = Y] 73 sin 2 (7r?7n). 

n=l M 



(19a) 
(19b) 

(20a) 
(20b) 

(21) 



Averaging the Gibbs-Thomson condition over individ- 
ual lamellae, we obtain the mean undercooling of the 
solid-liquid interface: 



(AT) a = m a ((c) a - c E ) + rha (S) a + T a (K) a (22a) 
(AT) = -m p {{c) fi - c E ) + rh (c) p + Tp (K) p , (22b) 

where the average values (AT),,, (c) s , and (K) s (s = 
a, f3) are defined by expressions analogous to Eqs. ( |l9| ) 
for (c) s . The averages for the composition, (c) g , and for 
the curvatures (K) s , are 



2AAc , . 
\c) a - c E = Coo + -Bo H Jr~~ 

(c) ^ - c E = Coo + B - — —P{r)) 



Ki -vY 



(K) a =^sin^ 
(K) g = 



P (1 - T))X 



sm ( 



(23a) 
(23b) 

(24a) 
(24b) 



The last step is to apply the condition that neighboring 
lamellae should grow at equal undercooling, (AT) = 
(AT) j3 . This determines the only degree of freedom left 
in the problem: the eutectic boundary layer Bo. The 
solution is 



B = -(coo - c E ) + 
2XP(r])Ac 

- 1 ( O^L 



rn - m C) 



ks m 



1-7? 



m 



1 



V 



Ac V 



kn k a 



(25) 



There are two terms that are not present in the binary 
eutectic. The second term on the RHS of Eq. ( |2^ ) gives, 
according to Eq. (|J), the shift of the eutectic composi- 
tion corresponding to an impurity concentration CoojkE- 
It is remarkable that, even for a system that started from 
an initial state with three-phase equilibrium, a eutectic 
boundary layer must develop in order that the composi- 
tion condition for three-phase equilibrium at the trij unc- 
tion points can be met. The second new term is the last 
term in brackets on the RHS of Eq. (p5f), involving the 
difference of the partition coefficients. This term is due 
to the unequal rejection of impurities into the liquid. 

The interfacial undercooling as a function of the lamel- 
lar spacing can be written in a form very similar to the 
JH result for a simple binary eutectic: 



AT = M^ 



1 



-AT 



A 



•^min 



(26) 



The first term, according to Eq. (m , gives the undercool- 
ing of the point in the eutectic valley corresponding to 
an impurity concentration Coo/kE- The minimum under- 
cooling AT m i n and the corresponding spacing A m i n are 



AT ■ — 



4Ac 



ry(l - -q) \m a + mp J 



(27) 
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with 



(1 - i])T a sm9 a r)Tp sin Op 



P(f?,Coo) = 



m a Ac 



/ 



At 



mpAc 



(k/3 - kc^j 



in a mp 

(l-r? 7? 

m a mp 



(28) 



(29) 



(30) 



Note that, as in a binary eutectic, we have Ar m ; n <• 

and A m in ~ In the special case k a — kp, we re- 

cover the classic JH result. The two A-dependent terms 
on the RHS of Eq. ( |2^ ) represent the effects of diffusion 
and surface tension. For finer lamellae, the diffusion be- 
tween adjacent lamellae is faster, and the undercooling 
due to the concentration term in the Gibbs-Thomson re- 
lation is smaller. This gives the term proportional to A. 
On the other hand, for finer lamellae the average curva- 
ture is higher due to the constraints at the trijunction 
points, leading to the term proportional to A -1 . 

For our subsequent stability analysis, we will now sim- 
plify the problem. The main effect of the impurities, re- 
jected by both solid phases, is the buildup of the impurity 
boundary layer of amplitude Bq. The diffusion of impu- 
rities between neighboring lamellae leads to corrections 
in A m j n and AT m j n . In the dilute limit, however, where 
Cqc -c Ac, these corrections are small. In addition, the 
partition coefficients for ternary impurities are often close 
to zero, and we have k a — kp <C 1. Thus, it seems well 
justified to make the approximation j?(t7, Coo) f=a P(rj)/l 
in Eq. (^TJ) and to drop the term involving the difference 
of the partition coefficients in Eq. (|2^) . This means that 
we neglect the interlamellar impurity diffusion modes, 
B n = for ri 7^ 0, which is equivalent to the assumption 
of equal impurity partition coefficients, k a = kp = &e- 
For two very different impurity partition coefficients, it 
might be necessary to go beyond this approximation and 
to include the interlamellar impurity diffusion. 

We will also assume equal impurity liquidus slopes rh a 
and mp for most of what follows. Then, the eutectic 
composition and the magnitude of the eutectic boundary 
layer, described by the coefficient Bq, do not depend on 
the impurity concentration. We will briefly comment on 
the general case at the end of Sec. VI. 



IV. STABILITY ANALYSIS 

The DL method is a perturbation analysis around the JH 
steady-state solution. The fundamental variables in this 
approach are the coordinates of the trijunction points, or 
more precisely the departure of these coordinates from 
their steady-state values. As illustrated in Fig. EL the 
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FIG. 4. Sketch of a perturbed lamellar interface, showing 
the displacements y?, yf, and 



coordinates of the trijunction points of the jth lamella 
pair are written as 



= j\ + yf(t) tf = (j + r,)\ + yf(t) 



(31) 



(32) 



We assume the system to consist of a total number of N 
lamella pairs, and use periodic boundary conditions. For 
convenience, we define a dimensionless eutectic concen- 
tration field by 



ce 



and a dimensionless impurity concentration field 



(33) 



(34) 



with Ac = Coall/kE — 1). 

Let us outline the strategy of the DL calculation. For a 
slowly evolving, slightly perturbed interface, the Gibbs- 
Thomson condition for local equilibrium remains satis- 
fied. The deformation of the front modifies the local 
curvatures and concentrations. The perturbed diffusion 
fields can be calculated using the mass conservation con- 
ditions. For small displacements, all resulting expressions 
are linearized in £| and y|. To close the set of equations, 
DL use Cahn's hypothesis, which gives additional rela- 
tions between the £'s and the y's. Inserting all these ex- 
pressions into the Gibbs-Thomson condition allows one 
to determine the stability spectrum of the interface. 

The new element here is the ternary impurity diffusion 
field. To obtain the modifications it generates in the 
stability spectrum, we have to introduce the procedure in 
more detail. As in the JH calculation, a complete solution 
of the free boundary problem is out of reach, and we 
use quantities that are averaged over individual lamellae. 
The average undercooling of a lamella is written as 



AT? (t) = AT + ST?(t) (s = a, /3), 



(35) 
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with AT, the steady-state value, given by Eq. (p6|). To 
express the 5T?(£)'s, we have to calculate the deviations 
of curvature and concentrations from their steady-state 
values, that is: 

S (K); (t) = K[C(x, t)]dx - (K) a , (36a) 

X 3 ~ X Jx 7 

S (10? (*) = - 1 - K[«x, t)]dx - {K) p , 

X 3+l ~ X J Jx i 

(36b) 



1 f 3 

S (u)°? (t) = — / 3 u(x, ((x, t),t)dx - (u) a , 

x 3 - ■>' 

(37a) 



be described in more detail below [see Eq. Fol- 
lowing this step, Eqs. (^) become a system of 2N linear 
equations for AN variables and their time derivatives. To 
complete this system, we must specify how the trij unc- 
tion points react to deformations of the growth front. 
Following DL, we will use Cahn's hypothesis and assume 
that the trijunctions always grow perpendicular to the 
eutectic solidification front, which yields the conditions 

# = -(£/-<?-l)VA W 

Vj = -{Zj+i ~ ij)v P /\ (40b) 

where the dot denotes the time derivative. 

We now transform the problem into an eigenvalue 
equation by analyzing it in terms of time-dependent 
Fourier modes. We write 

= X s k exp(ifcAj + ut) (s = a, 0) (41a) 



B 1 f Xj+1 
8 ( u )j (*) = — 8 / „ u ( x > t), t)dx - (u)g , 

X ?+l ~ X 3 J *i 

(37b) 

with the equivalent quantities for the ternary impurity 
field being obtained by replacing u by u in the last two 
equations above. The steady-state values are obtained 
from Eqs. (|2C| ) and ( |23| ) using the changes of variables 
( |33| ) and (|34|). Furthermore, we define the average z po- 
sition of a lamella by 

<0J=~(^+^) 08a) 

(Of = i(^+i + «/)- (38b) 

In terms of these local averages, the linearized Gibbs- 
Thomson conditions read: 

5T?(t)= -G(t)°(t) 

= m a S (u)°? (t) + rh a S (u)° (t) 

+ T a S (K)° (t) (39a) 

STf(t)= -G(0?(t) 

= - m 5 («)J (t) + (fi)J (t) 

+ i>* <#)?(*)• (3%) 

The next step is to express the averages in the two equa- 
tions above in terms of the displacements £| and y|. For 
the curvature terms, the procedure is straightforward, 
but for the diffusion fields one needs to introduce a piece- 
wise linear representation of the interface shape, as will 



y° =Y£exp{ik\j+cjt) (s = a,0), (41b) 

where the allowed wave vectors k are integer multiples of 
2-k/NX and lie in the interval [— tt/X,tt/X\. In the limit 
of an infinite number of lamellae, N — > oo, we recover a 
continuous spectrum. The growth constraints (f4(i| ) give 
then immediately 

uY a = _?^P e -ikx/2 sin(A . A/2 )xf (42a) 
A 

uY p = _2wp eikx/2 sin(fcA/2)x «_ (42b) 
A 

This allows us to eliminate the coefficients Y k a and Y k , 
and the only unknowns left in the problem are X k , X?, 
and u>. It is useful to write each of the terms appearing 
in the Gibbs-Thomson condition in the compact forms 





ikXj-\-ujt 


£ G s ' s> 


\k,u)X( 


(43) 




gik\j-\-ujt 






(44) 




ikXj-\-Ljt 


u s ' s ' 

s'=a,{3 


{k,U3)Xi 


(45) 






^ U s ' s ' 


{k,u)xi. 


(46) 



Then, the conditions Q3S|) can be written as an eigenvalue 
equation: 

£ A' - A/ 0. (47) 

s'— a.,0 

where the matrix A is given by 
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G 



r Q o 
o r> 

rh a Ac 



K 



m„Ac 




-mpAc 



U 



rhpAc 



U. 



(48) 



The first three terms on the RHS of Eq. ( |4q ) are identical 
to those calculated by DL, whereas the last term is due 
to the presence of the ternary impurity. For a given wave 
number k, Eq. (E^) is fulfilled only for special values of 
lu. The solvability condition, 



dct A = 0, 



(49) 



gives the dispersion relations uj(k). 

The core of the problem is the calculation of the ma- 
trices appearing in Eq. (f48|). The first two, G and K, 
are relatively easy (see Ref. 0] and appendix A). The 
most difficult is the matrix U. We start by writing the 
perturbed diffusion field as 



u(x, z, t) = uq(x, z) + 5u(x, z, t), 



(50) 



where uq(x, z) is the steady-state solution. Next, 5u is ex- 
panded in Fourier modes, this time with the periodicity of 
the whole system, i.e. with wave numbers p = 2Tm/NX: 



8u(x, z, t) 



p 



b p exp[ipx — q p (z — z) + u>t\. (51) 



The p's, unlike the fc's, range from — oo to oo, because 
the diffusion field is continuous. To simplify our calcula- 
tions, we use the quasistationary approximation for the 
diffusion equation, i.e. we drop the time derivative in 
Eq. (|7a|). Physically, this means that we assume that 
the diffusion field adjusts instantaneously to any change 
in the interface configuration. For a perturbation such as 
that given by Eq. (|5l"|), this is justified if \oj\ <C Dp 2 and 
\uj\ <C Dp/l, conditi ons t hat will be checked a posteriori 
at the end of Sec. VD. In the range of wavelength of 
interest here, i.e. for A < 2n/p < I, we find that these 
conditions are generally well satisfied. Within this ap- 
proximation, the constants q p are given by the analog of 
Eq. @. 

To obtain the Fourier coefficients b p , we proceed as in 
the steady-state calculation and insert the expansion in 
the mass conservation condition. To make the problem 
tractable, the actual interface shape £(x, t) is replaced by 
the piecewise constant function 



COM) 



(0- 
(0" 



Xj < X < Xj 

Xj <C x <C Xj^,-^. 



(52) 



One needs to be careful because the gradients of the 
steady state concentration field diverge at the trijunc- 
tion points. Details can be found in DL's article. Fi- 
nally, the result is inserted in Eqs. (|37]). There are two 
types of contributions to linear order in £| and yj: the 
steady-state diffusion field is evaluated at the position 



of the perturbed interface, and the perturbed diffusion 
field is taken at the location of the steady-state inter- 
face. These different terms, containing sums over the 
steady-state interlamellar diffusion modes, lead to quite 
complicated expressions, summarized in appendix A. 

The calculation of U is somewhat easier because, with 
the assumption of equal solute partition coefficients, the 
steady-state impurity diffusion field is translationally in- 
variant along x, and hence the horizontal displacements 
Uj drop out of the calculation. It would be possible to 
include the interlamellar impurity diffusion in the expres- 
sion for U by following the lines of the DL calculation for 
U. We expect, however, as argued in the preceding sec- 
tion, that this would only lead to minor corrections. 

The perturbed impurity diffusion field is written as 

u(x, z, t) = uo(z) + 8u(x, z, t) (53) 
and expanded in Fourier modes, 

5u(x, z, t) = b p exp[ipx — q p (z — z) + uit]. (54) 
p 

To calculate the constants q p , we use the quasistation- 
ary approximation of the impurity diffusion equation. In 
terms of the dimensionless impurity field {£, the continu- 
ity equation at the interface takes the form 



D 



di 



(v p + t) [(l-k E )u- 



(55) 



We want to keep only terms that are linear in the dis- 
placements £| or their time derivatives. Such terms come 
from several sources: the time derivative of the correc- 
tions Su in the diffusion field and its gradient, and from 
evaluating the steady-state diffusion field at the new in- 
terface position £(x,i). The equation of order 1 in the 
displacements becomes 



-D a2 "° 



dz 2 
C[(l 



z—z 

- k E )u Q (x, z) + k E ] 



v p (l - k E ) -t— 
oz 



(C-z) 



+ v p (l — k E ) Su(x, z, t). 



(56) 



We now insert the Fourier expansion (^4|) in the above 
equation, multiply both sides by exp(—ip'x)/NX and in- 
tegrate over x from to NX. With uq(x, z) = exp(— 2(z — 
z)/l) this leads to 



q p - j(l - k E ) 



NX 



4k E 1 
I 2 NX J0 

1J_ f NX 

+ D NX J 



- ipx {(x,t)dx 
- %px C,{x,t)dx. 



(57) 
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To perform the integrals over x, we use the piecewise con- 
stant expression Eq. for £(x, t). As (—z is already of 
order £, we can neglect the y's in the integration bound- 
aries. For £| given by a Fourier mode of wave number k 
according to Eqs. mw, we obtain 



pi] A 
p(l- 



4fc £ 
f 2 



X 
p 



-iprjX/2 



-ipX/2 <rf3 



(58) 



where 



X„ 6 



JV-l 



iV ^ 



-ipXj—ut cs 
>7 



= /J ^>,fc+27rn/A-^jb ( s = /?)• 



(59) 



Note that, even if we start with a set of £'s given by a 
single Fourier mode of wave number k, the use of the 
piecewise constant interface shape of Eq. ([32]) induces 
perturbations in the diffusion held at all wave numbers 
shifted with respect to k by a multiple of 2tt/X. This 
effect is unavoidable if we want to replace the full free 
boundary problem by equations for a discrete set of vari- 
ables. 

The next step is the calculation of the average concen- 
trations in front of each lamella: 



TjX 



jA 

(&>?(*) = 



\ dz 



+ 5u(x,z,t)) dx (60a) 



1 



(Oi + Su(x,z,t) ) dx. (60b) 



(l-r?)A 

'0"+'?)A V dz 

The hrst term inside each integral is due to the displace- 
ment of the interface in the steady state diffusion field 
and the second arises from the perturbed diffusion field 
evaluated at the steady-state interface position. The final 
result for the matrix U is 



fjfl,f> 



1 1 2k E X luXI \ 
1 7]! \ I 2D J 

x &(k,»7) + §2{K,T)) 
1 1 / 2k E X uXf 



I (1-T))l\ I 



2D 



Si(k, 1 - rj) + SZ(k, 1 - ??) 



(61a) 
(61b) 



(61c) 
(61d) 



where U a > a * denotes the expression obtained from U Q,a 
by complex conjugation of all quantities except for u> (for 
u> real, this is the usual complex conjugation), and we 
have defined 



K = kX/2iT 



r = 1/1 = D/D 



(62) 



(63) 



Pn{K) 



r 2 Pe^ 



47r 2 (n 



n— — 00 



sin 2 [7r?7(n + «)] 
7r 2 (n + k) 2 p„(k) 



rPe + 2rPek E (64) 



(65) 



-in[n-\-K) 



sm[iri](n + k)} sin[7r(l -rj){n + k)} 

These notations have been chosen in analogy to some of 
DL's results (see appendix A). The ratio r of the eutec- 
tic and impurity diffusion lengths is usually close to one, 
and k is the dimensionless wave number. For small Peclet 
numbers, and perturbation wavelengths much larger than 
the lamellar spacing (k <C 1), the sums S\ and S% are 
dominated by the term with n = 0. In this limit, we 
can neglect all the other terms in the sums, which cor- 
responds to keeping only the Fourier coefficient b p with 
p = k, and hence to a single- mode approximation of the 
perturbed impurity diffusion field. When k is larger, 
however, and in particular near the "Brillouin zone" 
boundary, k = 0.5, we have to consider the full sums. 
To obtain the stability spectra, we must now combine 
this result with DL's calculations for the other matrices 
and solve the characteristic equation for u>. 



V. SYMMETRIC PHASE DIAGRAM AT 
EUTECTIC COMPOSITION 

A. Stability spectrum 

In general, the characteristic equation of the stability 
spectrum is a polynomial of degree four in to with real co- 
efficients. The solutions could be obtained algebraically, 
but this method leads to complicated expressions that are 
difficult to interpret. For this reason, we will restrict our 
attention in this section to the special case of an alloy of 
eutectic composition in a model system where the phase 
diagram is symmetrical about the eutectic composition. 
In this case, the characteristic equation can be factored 
into two quadratic equations, thus greatly simplifying the 
interpretation. As we shall see when we treat the general 
case in Sec. VI, all qualitative features of the instability 
are already contained in this special case. 
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There are general relations between the elements of the 
matrix A, due to the existence of two planes of mirror 
symmetry in the steady state, one in the middle of each 
type of lamella (a or (3). Hence we can change the sign 
of k and relabel the trijunction points without affecting 
the final result. This leads to the relations 



A a,/S = A av 



A /3,/3* 



(67a) 



(67b) 



Here, an asterisk again denotes complex conjugation of 
all quantities except lu. If we consider a model cutcc- 
tic with a completely symmetric phase diagram, i. e. 
m a = mp = m, rh a — rhp = rh, T a = Tp = T, 
9 a = Op = 9, and up = —u a = 1/2, at its eutectic 
composition, i] = 1/2, we have in addition A a ' a = AP'^ '. 
The solvability condition, detA = 0, can then be fac- 
tored into two equations, 



Re (V jfeA/4 A Q < Q ) = and Im ( e - ihX / i A a > a 



0. 



(68) 

both of them quadratic in u>. 

To proceed, we will rewrite the equations in a dimen- 
sionless form. For the sake of subsequent generalization, 
we will give expressions for the parameters that are valid 
for any phase diagram. We define 



M = 


m a mp 


(69) 


(m a + mp)/2 


n = 




(70) 




Gl 2DG 


(71) 


9 = 


MAc VpMAc 


w 


toAc 
MAc 


(72) 


A 




(73) 



Here, M is a mean liquidus slope (for the symmetric 
phase diagram, M = to); 1/Ref2 is the distance along 
the z direction that the interface needs to travel (in 
units of lamellar spacing) for the amplitude of the per- 
turbation to grow by a factor of e, and 27r/Imf2 is the 
length traveled by the interface during one oscillation 
cycle. The parameter w is the ratio of the two freezing 
ranges and is proportional to the impurity concentration 
since Ac = (1/fc^ — l)coo. 

It may be useful to comment on typical experimental 
values of the dimcnsionless parameters. The best studied 
system in this context is the organic eutectic CBr4-C2Cl6 
used originally by Hunt and Jackson [Q. This system 
contains naturally residual gas that acts as an impurity, 
leading to colony formation. For this system, Akamatsu 
and Faivre [Q have estimated w to be of order 0.1, with 
a distribution coefficient k E ~ 0.02. The Pcclct number 
Pe is typically between 0.01 and 0.1 in the low-velocity 
regime used to investigate colony formation. 



After multiplication by 1^1/ (MAc), the first of Eqs. 
(|68| ) becomes, expressed in the above dimensionless quan- 
tities (see appendix A for more details), 

„ q TTK 

= fl- cos — 
2 2 

2P(n) 

H i^- sin(7TK/2) sin(TTK) (Vlcote - 2/rj) 

7]A Z 

+ fffic(V" rK / 2 J7 2 Q (K,77) 
+ 2sin(7TK)Re f ie t7TK/2 U 3 (K, rj) 



(74) 



„ 7tk wrQ , 

wrfl cos — H (2k E rPe + CI) 

2 rj 

§j_ ( K) n) cos ^ + Re (V" rK/2 S , *( K , ??)) 



where we have chosen to display rj, for clarity, although 
the factorization is only possible for r/ = 1/2. The first 
term of Eq. (|74| ) arises from the matrix G; the factor 
cos(7tk/2) is simply due to the averaging over the two 
trijunction points limiting a lamella. The matrix K con- 
tributes the next two terms, proportional to A~ 2 . The 
first, containing f2, describes the change of curvature due 
to the bending of the interface over a large scale, and 
is therefore equivalent to the capillary term in the di- 
lute binary alloy problem. The second gives the change 
in average curvature upon variation of the local lamel- 
lar spacing. All terms containing the functions U"(r], k), 
defined in appendix A, are due to the eutectic diffusion 
field. Finally, the terms proportional to w arise from the 
matrix U. 

Equation (|74] ) is exact and can be solved numerically. 
But it is also useful to simplify this equation in order to 
render the physical interpretation of the instability more 
transparent. To this end, let us group the terms with 
equal powers of fl and rewrite Eq. ( 74 ) as 



a(n)Q 2 - b{n)n + c{k) = 0. 



(75) 



To obtain a simplified expression, we expand the coef- 
ficients a(ft), b(n), and c(n) in powers of k. Details on 
this procedure can be found in appendix B. It turns out 
that an expansion up to order n 2 is sufficient to obtain a 
satisfying agreement with the direct numerical solution 
of Eq. f74|). We obtain 



Po( K ) 

27r 2 P(l/2)cot0 
N 2 



a(n) = 2P(l/2) - 

6(k) = wr — g/2 

2wPer 2 kE 

c(«) = 8tt 2 P(1/2) (l - « 2 , 



(76a) 



i?0 

~2~ 



K 2 



(76b) 
(76c) 
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FIG. 5. Stability spectra of the symmetric eutectic at 
A = 1, Pe = 0.01, g = 1, r = 1, k E = 0.05, and 9 = 45°, 
(a) without impurity (w — 0), (b) with a small amount of 
ternary impurity (w — 0.01). Full lines: real modes, dashed 
lines: complex modes. 



where Po(k) is defin ed by Eq 
given by Eq 



(Bll 



(p4|), and Rq is a constant 
We note that, to obtain these ex- 



pressions, we have only kept the leading order terms of 
each of the contributions in Eq. (jTJ) and dropped several 
terms of order |re| and re 2 that turn out to give negligible 
contributions at the onset of instability for the reasons 
detailed in appendix B. The simplified stability spectrum 
defined by the equations above will be used below to de- 
rive simple analytical expressions for the onset velocity 
and wavelength. Moreover, it will allow us to identify the 
terms that contribute to the effective surface tension in 
the long- wavelength free boundary formulation presented 
in Sec. VI. In the rest of this paper, the results based on 
this simplified spectrum will be systematically checked 
against the direct numerical solution of Eq. ([74|) . 



B. Binary versus ternary eutectic: the limit re 







Let us first examine the limit re — > of the stability 
spectrum, which governs the relaxation to steady-state 
growth after an infinitesimal translation of the solidifica- 
tion front along the z axis. Although this is not the limit 
of interest for morphological instability, it is worth a brief 
discussion to understand more fully the subtle effects of 
the ternary impurity on the complete stability spectrum. 

Since, as mentioned above, the complete eigenvalue 
equation is of degree four in O with real coefficients, there 
are four branches of the dispersion relation, and modes 



can be real or occur in complex conjugate pairs. In the 
completely symmetric case, each of the two equations 
( |68| ) gives a pair of branches. To calculate the stability 
spectra, we must first evaluate numerically the sums oc- 
curring in the functions U n (n, rj) to obtain the coefficients 
of the polynomial in Q, which can then be solved for each 
re. Fig. H shows a comparison between the stability spec- 
trum of a binary eutectic [w = 0) and a eutectic with 
a small amount of impurities, "small" meaning that we 
stay far below the threshold of instability. We have cho- 
sen the parameters g and A in order to reproduce, for the 
binary eutectic, Fig. 4(b) of DL's article. Without im- 
purities, we can distinguish two types of branches. There 
is always a characteristic diffusive branch, which is real 
and satisfies ft oc 2?(A)re 2 for re — > 0: the spacing can 
be locally adjusted by "diffusion" in A space [Q. This 
branch is related to the long-wavelength lamella elimi- 
nation instability for A < A m j n : the effective diffusion 
coefficient T>(A) is negative for A < 1. The companion 
of this branch is also real, with strongly negative growth 
rates. The second pair of branches is complex for small 
k. This mode gives rise to the 2A-0 instability for suffi- 
ciently off-eutectic compositions: it becomes complex up 
to the "Brillouin zone" boundary k = 0.5 and its real 
part becomes positive. 

When we add a small amount of impurities (w = 0.01), 
we find that the diffusive and oscillatory branches are 
nearly unaffected; however, the companion of the dif- 
fusive branch undergoes a drastic change. Physically 
speaking, this strong reaction to a seemingly small per- 
turbation is due to the introduction of a new conservation 
law. In the binary eutectic, this mode describes the re- 
laxation of an interface by the motion of the trijunction 
points with respect to the temperature gradient (in the 
z direction). At the eutectic composition, this relaxation 
involves only short range interlamellar diffusion and is 
therefore fast. On the other hand, for a flat interface, the 
impurities must diffuse over a length scale of the order of 
the diffusion length to achieve relaxation. This leads to 
much slower decay rates for small wave numbers. 

To see more formally how this change arises, let us 
consider Eqs. ( fT6| ) in the limit re — ► 0. Without the 
impurity terms (w = 0), a(n) and 6(re) stay finite in the 
limit re — ► 0, whereas c(re) oc re 2 . Hence fl oc re 2 leads to 
a balance between the last two terms of Eq. ( |75| ) . In the 
impurity terms, we have po(n) — > 2/c^rPe for re — > 0, and 
for re = we obtain 



2P(l/2) 



2k E Pe 



fi 2 -fft = 0, 



(77) 



yielding the solutions, fi = 0, for the diffusive branch 
and, SI = —g/[AP(l/2) + w/(kEPe)], for its compan- 
ion. In the latter expression, the impurity contribution 
is dominant for small Pcclct number and small impurity 
partition coefficient, and hence this branch is strongly in- 
fluenced by the addition of impurities. As the terms pro- 
portional to po( K ) _1 also appear in the classical Mullins- 
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Sekerka analysis of a monophase solidification front, we 
will hereafter refer to this branch as the MS branch. 

The oscillatory branch is little affected by the addi- 
tion of impurities. The reason is that the second of Eqs. 
( |68| ) does not contain the MS terms at n = 0. We will 
discuss the relation between long- and short-wavelength 
instabilities in more detail in Sec. VII. 

The derivation of Eqs. ((74|) and (|7^) is based on the 
quasistationary approximation of the two diffusion equa- 
tions. This approximation relies on the assumption that 
the wavelength of the perturbation is smaller than the 
diffusion length, and hence breaks down for k -C Pe. In 
the framework of the DL formalism, however, the cal- 
culation becomes extremely tedious if this assumption 
is relaxed and the results of this calculation will not be 
displayed here. The effective interface formulation to be 
presented in Sec. VI, however, easily allows to relax this 
assumption and to include the dynamics of the diffusive 
boundary layer. As a result, the MS branch of the spec- 
trum becomes complex for k <C Pe, corresponding to 
an oscillatory relaxation of the interface to steady-state 
growth driven by oscillations of the impurity boundary 
layer that is already well-known for a monophase front. 
As we shall see below, the morphological instability lead- 
ing to colony formation involves only modes with K>Pe 
for which the quasistationary approximation is valid. 



C. Onset of instability 

Let us now examine the onset of instability and compare 
our findings to the well-known results for dilute binary 
alloys. For the one-sided model, the constitutional super- 
cooling (CS) criterion is fairly accurate. This criterion 
states that a monophase solidification front is unstable if 
the diffusion length is less than twice the thermal length 
It = toAc/G. In our dimcnsionlcss variables, this is 
equivalent to g < 2wr. The Mullins-Sekerka analysis 
shows that the actual critical velocity differs from 
CS by corrective factors that are usually small. The crit- 
ical wavelength A c at the onset of instability scales as 
A c ~ (dolrl) 1 ^ 3 , where do is the capillary length. 

Let us briefly comment on some consequences of this 
scaling. From Eq. ( p8f ) we can deduce that the eutectic 
spacing Amin scales as Amin 

(cIqI) 1 / 2 . Therefore, for low 
Peclet numbers we always expect A m ; n <C A c , and hence 
we can consider the limit of small k for the determination 
of the onset of instability. On the other hand, A c will 
always be smaller than the diffusion length, and we have 
M/2tt ^> 1, or equivalently k ^S> Pe. Therefore, we may 
use the simplification Po{k) rs 2~k\k\ in Eqs. (76). 

The occurrence of unstable modes is determined by the 
behavior of if wr — g/2 is positive and large enough, 
becomes positive for a certain range in k. We want to 
determine the critical value g c of the parameter g where 
the first unstable mode occurs, and the wave number k c 
of this mode. The two solutions of the quadratic equation 



|) are 



b(n) ± ^b 2 ( K ) -4cl(k)c(k) 
2o(k) ' 



(78) 



As the product a(/t)c(/t) is always positive and finite, the 
solutions become complex when &(k) tends to 0, and we 
have Refi = &(«)/2a(«). Hence the two conditions to ob- 
tain g c and n c are simply b(n) = and db[K)/dn — 0. Let 
us rewrite the first of these conditions in the dimensional 
quantities and divide through by wr. The result is 



1 



GD VI cos 9 k 2 i? A 2 k 2 



rhAcv„ 



mAc 2 ATT 2 wr 2 



2k E 

ki 



0. (79) 



There are two terms proportional to k 2 . The first is the 
surface tension term which stabilizes the interface. It is 
analogous to the surface tension term in the monophase 
MS spectrum, except that it is multiplied by a geometri- 
cal factor cos 9. This factor is present because the eutec- 
tic interface is made of an array of arcs, each one linking 
two trij unctions, which renders the front less stiff than 
a flat monophase interface. The second term arises from 
the eutectic diffusion field. It can be directly interpreted 
by noting that the above equation becomes identical to 
the result for a monophase planar interface if we define 
the effective capillary length 



d 



Tcosfl i? n PeA 



toAc 



An 2 w 



(80) 



The above expression implies that the interlamellar eu- 
tectic diffusion field has a stabilizing effect. This is rather 
remarkable since it implies that the two diffusion fields 
(associated with the eutectic components and the ternary 
impurities, respectively) play antagonistic roles in the in- 
stability. 

From the condition db(n)/dhi — we obtain the ex- 
pression 



kt = 



2k E 
d l 2 



(81) 



for the critical wave number. Furthermore, substituting 
this expression in Eq. (|79|), we find 



GD 



1 



toAc 



1 



(do 



kl/2l 



A 1/3 



(82) 



for the critical pulling speed. This expression is equiva- 
lent to the constitutional supercooling criterion (ves = 
GD /to Ac) up to the capillary correction in the denomi- 
nator of the second term on the RHS of Eq. ( p2| ) . This 
correction is typically negligibly small, except for very 
low impurity concentrations where the effective capillary 
length becomes large. Consequently, as for a monophase 
front, there is a critical impurity concentration below 
which the eutectic front is completely stable. 
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FIG. 6. Stability spectra for varying g and A and w — 0.1, 
Pe = 0.01, r = l,k E = 0.05 and 9 = 45°. (a) g = 0.15, A = 1; 
(b) g = 0.05, A = 1; (c) g = 0.15, A = 1.1; (d) g = 0.05, 
A = 1.1; (e) g = 0.15, A = 1.5; (f) g = 0.05, A = 1.5. Solid 
lines: real modes; dashed lines: complex modes. Note that 
we show only the diffusive and MS-branches in the region of 
small k where they are unstable. The inset in (a) shows in 
more detail the part of the spectrum at small k. 



The above results show that the instability of the eu- 
tectic interface is qualitatively similar to the standard MS 
instability of a monophase front, as far as the expressions 
for fc c and v c are concerned, up to a renormalization of 
the surface tension as described by Eq. (80). One main 
difference, however, is that the branch of the spectrum 
that becomes unstable is complex, which suggests the 
existence of oscillatory patterns with a wavelength much 
larger than the lamellar spacing. To decide whether such 
patterns are observable experimentally, we need to ex- 
amine next how this complex branch evolves above the 
onset of instability. 



D. Above onset 

For v > v c , the solutions for £1 still satisfy the quadratic 
equation (|75|). The nature of the mode, whether complex 
or real, is determined by the sign of the discriminant, 
b 2 {n) — 4a(K)c(ft). The most important parameter, be- 
sides g, that controls the discriminant is A, because c(k) 



strongly depends on A near A = 1. We are interested 
only in the case A > 1, since spacings below A m ; n are in- 
trinsically unstable. For A near 1, c(k) is small, and real 
modes should appear. In contrast, for larger A, complex 
modes should be favored. 

To check this prediction, we plot in Fig. || a series of 
stability spectra, calculated using Eq. (74), for varying 
g and A. We chose w = 0.1, which, using the constitu- 
tional supercooling criterion, gives a critical temperature 
gradient of gcs = 0.2. We display spectra for two values 
of 5, one close to onset (g = 0.15), and one far above the 
onset (g — 0.05). 

Let us first comment on the structure of the spectrum 
at small values of n. The inset of Fig. ||(a) shows that the 
growth rate still satisfies Eq. ( |77| ) for k = 0, but now as 
k grows the MS branch and the diffusive branch meet to 
form a complex conjugate pair, because b(n) approaches 
zero and the discriminant becomes negative. When k 
grows further, the real part of Q becomes positive, and 
the branch may stay complex or split again in two real 
branches, both with positive growth rates. For A = 1 
[Fig. ^ (a) and (b)], the spectra always exhibit real modes 
above the onset, and the mode with the maximum growth 
rate is real for both values of g. For A = 1.1 [Fig. ^ 
(c) and (d)], near the onset only a narrow band of real 
modes is present, and the fastest growing mode is com- 
plex. This changes at lower g: the k range of real modes 
has increased, and the fastest growing mode is real. Note 
that in both spectra there are two maxima of the growth 
rate, one real and the other complex. At a certain value 
of g, the two maxima are of equal height; at this point, 
we have a finite jump in the wave number of the fastest 
growing mode when g is varied, and two perturbations 
of different wave numbers grow with the same rate. Fi- 
nally, for A = 1.5 [Fig. |TJ (e) and (f)], there is a range of 
g where the spectrum is entirely complex, and only far 
above the onset real modes appear. The fastest growing 
mode is always complex. 

We can summarize these results in a diagram that 
shows the nature of the spectrum in the plane (A, g) (Fig. 
0) . For this diagram, we normalize the temperature gra- 
dient by its critical value according to the constitutional 
supercooling criterion, gcs = 2rw. Note that we have 
9 /9CS = vcs/vp, where vqs = GD/rhAc is the critical 
velocity. We classify the spectra into four categories, ac- 
cording to the occurrence of maxima. We denote by R 
spectra with only a real maximum (Fig. ^|a); spectra 
with a real and a complex maximum are denoted by RC 
when the fastest growing mode is real (Fig. |^d), and CR 
when it is complex (Fig. ^3). Finally, the entirely com- 
plex spectra are denoted by C. Figure [?] shows for which 
parameters we can expect predominantly real or com- 
plex modes; around the line between RC and CR, we can 
have the competition of two different modes. This dia- 
gram was determined using Eq. ( |7^ ) with the parameters 
w = 0.1, Pc = 0.01, and ks — 0.05. Increasing the con- 
centration of impurities (increasing w) favors real modes: 
all curves are shifted to the right and to the bottom in 
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FIG. 7. Structure of the dispersion relation as function of 
the parameters A = A/A m i n and g/gcs = vcs/v p . The onset 
of instability (constitutional supercooling) is at the bottom of 
the diagram. The meanings of R, RC, CR, and C are defined 
in the text. 



the diagram. The Peclet number and the partition co- 
efficient influence the diagram only in the part near the 
onset of instability. 

To describe a spectrum in more detail, we may use 
several characteristic wave numbers: the wave numbers 
of the fastest growing modes, the limit between real and 
complex modes, and the wave number of the marginally 
stable mode n m . Figure ^ shows these quantities as a 
function of g for A = 1.1. As can be seen from Figs. ||c 
and |^d, for this value of A real modes become dominant 
only well beyond the CS threshold. 

To investigate systematically these quantities as func- 
tions of the control parameters, it is cumbersome to use 
Eq. ( |74| ) because of the sums involved in the functions 
U n . Therefore, it is convenient to obtain approximate 
expressions for the wave numbers of the marginal mode 
and the complex maximum using the approximate spec- 
trum defined by Eq. ( f75|) and Eqs. ([76]). We will show 
in Sec. VI that this simplified spectrum can be recovered 
from an effective interface approach, which applies to an 
arbitrary phase diagram and composition. It is therefore 
worthwhile to investigate the quality of this approxima- 
tion in the present symmetric case by quantitative com- 
parison with the exact spectrum of Eq. (Q). Figure ^ 
shows that, for the example of the last spectrum in Fig. 
||, the qualitative aspect of the spectrum is well repro- 
duced. We checked several cases and always found that 
the approximation shifts the marginally stable mode and 
the maxima to larger k. The error, however, never ex- 
ceeded about 30%. The imaginary part of f2 is very well 
approximated over the whole range of k. We also recal- 
culated the diagrams of Figs. and||, and found that the 
lines are slightly shifted, but the qualitative structure of 
the diagrams stays the same. It seems therefore valid to 
use Eqs. (|7^) for a general analysis. 

The marginally stable mode is always complex, and 
fcyn can be determined by the condition &(/c m ) = 0. In 




FIG. 8. Wave numbers of some characteristic points in 
the stability spectrum as a function of g/gcs for W = 0.1, 
Pe = 0.01, A;_b = 0.05, and A = 1.1. Solid lines with circles: 
real maximum, with crosses: limit between real and com- 
plex modes, with triangles: complex maximum, and with dia- 
monds: marginally stable mode. The dashed and dash-dotted 
lines are the approximations for the complex maximum and 
the marginal mode given by Eqs. (Ei]) and (§51), respectively. 



addit ion, far enough above the onset the term in 1/k in 
Eq. (76b) can be neglected. We find 



1 / 9\ 
— [wi — - 

a \ 2/ 



with 



a 



27r 2 P(l/2)cot( 

A 5 



Ro 



do, 



(83) 



(84) 



where Ro, as before, is given by Eq. (Bll). We can 
also give an approximate expression for the wave number 
of the complex maximum. For complex modes, Re fj = 
6(k)/2o(k). If we neglect the first term in Eq ( 76a ) for 
o(k), we obtain 



(85) 



The resulting curves for K m and n rnax are shown in Fig. 

For the real maximum, all terms in the quadratic equa- 
tion for fl have to be retained, and no simple expression 
for the wave number of the maximum can be obtained. 

Some remarks seem to be in order here to clarify the 
meaning of Eq. ( |83| ) in the dimensional variables. Since 
9/9CS = v cs/v p , the limit of high pulling speeds cor- 
responds to g — ► 0. In this limit, the structure of the 
spectrum varies very little with the pulling speed, be- 
cause both the characteristic length of the MS problem, 
\J dol, and the spacing A m i n vary as 1/ far above the 
onset. Hence, if we increase the pulling speed and at the 
same time change the lamellar spacing such that A stays 
constant, the only parameter that changes (except g) is 
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FIG. 9. Comparison of stability spectra obtained by the 
full calculation, Eq. (Q, and by Eq. (||) using Eqs. for 
the coefficients a, fe, and c. Symbols: full calculation (circles: 
real modes, crosses: complex modes). Lines: approximation 
(solid lines: real modes, dashed lines: complex modes). The 
parameters are the same as for Fig. ^(f). 



the Peclet number. But Pe appears in the problem only 
in the last term of Eq. (76b), which we have neglected in 
order to obtain Eq. (|83|). This term becomes important 
only at very high pulling speeds and leads to the absolute 
stability of the interface, as in the case of a dilute binary 
alloy. Substituting k = kX/2-K in Eq. j8^ ) shows that, for 
constant A, far above the CS threshold the wavelength of 
the fastest growing complex mode scales as \J dol, as for 
monophase solidification. On the other hand, in exper- 
iments on eutectics unstable states are usually reached 
from stable states by a sudden increase of the pulling ve- 
locity |2]],|28| . The lamellar spacing immediately after the 
jump is the same as before, but A m ; n and consequently 
A have changed. Hence, starting from the same initial 
state and varying the final pulling speed corresponds to 
a variation of both g and A. Equation (83) still applies, 
but in view of Eq. (^4|) no simple scaling with v p is ex- 
pected. A last remark concerns the dependence of n m on 
the impurity content. Taking the limit g — > in Eq. ( p^ ) 
shows that n m ~ y/w ~ a/coo ■ The reason for this behav- 
ior is that the effective capillary length, Eq. (|bo|), scales 
as do ~ w , whereas A m i n is independent of w. This 
means that the MS instability length increases with de- 
creasing w, whereas the lamellar spacing stays constant, 
and hence we expect the wavelength of the primary in- 
stability of a eutectic front to decrease with increasing 
impurity content. 



Finally, let us verify that the quasistationary approx- 
imation we have used to obtain the above results is jus- 
tified. To this end, we estimate the order of magnitude 
of the three terms in the diffusion equation Eq. (Q). 
For a perturbation Su of wave number k > l/l (re > Pc) 
growing at rate w, we have dtSu ~ coSu, d z Su ~ k5u, 
and V 2 5u ~ k 2 Su. In terms of the dimensionless vari- 
ables, the magnitude of these terms is |fi|Pe, rePe, and 
k 2 . From Fig. || we see that, for k > Pe, |0| < Bk 
with some number B of order unity, and hence the omis- 
sion of the time derivative from the diffusion equation is 
well justified for the range k > Pe of interest here. It is 
also possible to relax the quasistationary approximation. 
Then, the growth rate oj appears in the denominators of 
all the sums S n (n) and 5„(k). An analytic treatment be- 
comes impossible, but the equation for can be iterated 
numerically. We have checked that for the present range 
of parameters, the use of this complete calculation leads 
only to insignificant changes in the spectra. Note, how- 
ever, that for higher Peclet numbers it may be necessary 
to include this effect. 



E. Oscillatory modes 

The most interesting result of this analysis is evidently 
the existence of complex modes. To illustrate the type 
of microstructures these modes would generate, we have 
calculated the trajectories of the trijunction points for a 
particular example, using the definitions Eqs. ( |4l"| ) and 
the complex amplitudes X| and Y£ obtained from the 
eigenvalue equation. For the symmetric alloy at the eu- 
tectic composition, the symmetry between the a and /3 
phases gives immediately X^ — exp(ik\/2)X£ . The am- 
plitude and phase of X£ may be chosen arbitrarily, as 
this amounts to fixing the origins of the space and time 
axes. The Kf were then calculated using the growth con- 
straints Eqs. (|TJ). 

We chose the complex modes with a maximum growth 
rate of Fig. |(e) (to = 0.1, A = 1.5, g = 0.15). As for ev- 
ery point in the complex part of the spectrum, there are 
two "degenerate" modes with complex conjugate growth 
rates. The modes of our example have a wavelength of 
10 A (re = 0.1) and growth rates ft = 0.0347 ± 0.257i. 
One of these two modes is depicted in the left part of 
Fig. [l(]. We see that it has a "traveling wave" structure. 
The wavelength in the z direction, expressed in lamel- 
lar spacings, is 27r/Imf2, and the propagation velocity 
is v/vp = ImO/27TK. There are two complex conjugate 
modes: one propagates to the left, and the other to the 
right. As their growth rates are equal, we can create any 
superposition of the two, in particular a "standing wave" 
shown in the right part of Fig. [[o| 

The reason for the existence of these oscillatory modes 
is the interplay between the destabilizing impurity dif- 
fusion field and the dynamical response of the internal 
lamellar structure. A protrusion of the interface rejects 
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FIG. 10. Microstructures created by complex linear growth 
modes. Left: a single mode traveling to the left. Right: a 
standing wave mode obtained by the superposition of two 
complex conjugate modes. Growth direction from bottom to 
top. 



impurities more efficiently than a flat interface, and hence 
grows faster. But as the curvature of the front increases, 
the trajectories of the trijunctions are more and more 
curved, and the local spacing increases. This leads to 
a decreased efficiency of the interlamellar eutectic diffu- 
sion and, hence, the interface slows down. As a result, 
the protrusion grows back. The lamellar spacing, how- 
ever, still increases due to the geometric constraints, and 
the process overshoots, leading to a concave deformation 
of the front. This gives a geometric interpretation of the 
difference between traveling and standing waves: for the 
traveling wave, the lamellar spacings to the right and to 
the left of the protrusion are different, providing a driv- 
ing force for the propagation of the perturbation. For the 
standing wave, the interface shape and the spacing are 
"in phase" , and the perturbation oscillates without prop- 
agation. Hence, the resulting microstructure depends on 
the initial relation between interface shape and lamel- 
lar spacing. This implies that in an experiment, where 
the initial perturbations of interface position and lamel- 
lar spacing have no reason to be in a particular phase 
relation, one should observe all possible superpositions. 

There are two characteristic quantities (besides the 
wave number) related to an oscillatory mode: its prop- 
agation velocity, or equivalently the ratio of the wave- 
lengths in the z and x directions, 27TK/Imf2, and the ra- 
tio of its frequency and its amplification time, Imf2/Ref2, 
which determines how many oscillations should be ob- 
servable before the amplification leads to lamella elimi- 
nation. In Fig. [ll] we show both quantities as a function 
of g for A = 1.5, where complex modes always dominate. 
We used Eqs. (75) and (f76|), determined for each value 
of g the wave number of the fastest growing mode, and 
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FIG. 11. Ratio of oscillation frequency to amplification 
rate, Imf2/ReS7, and ratio of the wavelengths in the z and 
x direction, 27rft/ImS7, for the oscillatory modes with maxi- 
mum growth rates as a function of g/gcs- Other parameters 
are A = 1.5, w = 0.1, Pe = 0.01, k E = 0.05. 



calculated Refi and Imf2 at this point. Approximate val- 
ues could be obtained by using Eq. (|85l) to obtain Km.n.*- 
Figure 
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shows that the ratio of frequency and amplifi- 
cation rate diverges when we approach the onset. This is 
to be expected, as the growth rate vanishes whereas the 
imaginary part of Q remains finite. With decreasing <?, 
fewer and fewer oscillations are observable before lamella 
terminations occur. However, even at g = 0, more than 
two oscillations are completed before the amplitude of 
the mode has grown by a factor e, which means that 
such modes should be transiently observable. The am- 
plification ratio generally increases when w decreases or A 
increases. The ratio of the wavelengths is fairly constant 
and slightly increases when A decreases or w increases. 



VI. EFFECTIVE INTERFACE APPROACH 

The discrete analysis of Sec. IV can in principle be used 
to calculate the stability properties of a eutectic front 
for an arbitrary phase diagram and composition. As we 
have seen, however, the resulting eigenvalue equation is 
quite complicated. It is therefore advantageous to de- 
velop an alternate formulation of the stability problem 
by exploiting the fact that the instability wavelength is 
typically much larger than the lamellar spacing. The 
idea, therefore, is to consider the shape of the large-scale 
front instead of the actual lamellar interface, and to solve 
a modified free boundary problem for this "effective in- 
terface" , with boundary conditions that account for the 
effect of the underlying lamellar structure. 

It is useful to present this approach in two steps. In 
a first step, we write down the free boundary problem 
for the effective interface in the absence of surface ten- 
sion effects. This yields a rigorous long-wavelength limit 
where the expression for fl agrees up to order k 2 with 
the one obtained from taking the small n limit of the 
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full discrete spectrum. This expression also reduces, in 
the absence of a ternary impurity, to the one derived by 
Langer for a binary eutectic, and contains the long- 
wavelength instability leading to lamella termination for 
A < 1. In a second step, we introduce phenomenologi- 
cally the effect of surface tension guided by the insights 
of Sec. V. In the following, we will allow the volume frac- 
tion j] and the eutectic liquidus slopes m a and mp to be 
arbitrary, and we will only require for brevity of notation 
that rh a = ihp. We will briefly comment on the general 
case, where rh a ^ rhp, at the end of this section. 



A. Long wavelength limit 

We start by defining the effective interface as the con- 
tinuous curve £,(x,t) that interpolates between the dis- 
placements of the trijunction points and £f and a 
continuous field y(x,t) for the displacements along the 
x direction. To obtain the free boundary problem that 
governs the large-scale motion of this interface, we start 
by writing the diffusion equation for the ternary impurity 
in the liquid phase and the associated mass conservation 
condition at the phase boundary, which yields the equa- 
tions 



d t c= DV 2 c, 
-Dd n c = (v p + £)(!- k E )c, 



(86) 
(87) 



where k E is given by Eq. (|16|) and we have used in Eq. 
( |S7j ) the expression for the normal interface velocity in 
the moving frame, v n — v p + £, which is valid for small 
amplitude deformations of the interface. Next, we need 
a boundary condition for c on this interface. For this 
purpose, we note that lamellae can be assumed to grow 
locally in steady state as long as the interface deforma- 
tion is on a scale much larger than the lamellar spacing. 
Therefore, we can assume that for such deformations, c 
obeys locally the Gibbs-Thomson condition 



T = T F 



rhc- AT JH (X,v n ), 



(88) 



where the contribution of the eutectic structure to the 
interfacial undercooling is given by the Jackson- Hunt for- 
mula 
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(89) 



with AT min and A min given by Eqs. (|27j) and (gSJ) . Fi- 
nally, to complete the problem, we need to relate the 
local lamellar spacing A(a;, t) and the shape of the front 
£(x,i). This is done, as in Ref. |23]j , by noting that the 
local lamellar spacing is given by 



X{x,t) m A 1 + 



dy 
O.r 



(90) 



for y <C A where Ao is the unperturbed spacing. The field 
y can then be eliminated by using the geometrical con- 
straint that lamellae grow locally perpendicular to the 
solidification front, which, expressed in terms of the con- 
tinuous fields, takes the form 



dy{x,t) 
dt 



Up dx- 



(91) 



Equations (|86]-pl|) , together with the boundary condition 
c = Cqo far from the interface, define the free boundary 
problem for small amplitude and long-wavelength defor- 
mations of the effective interface. 

The stability spectrum can now be obtained by car- 
rying out a standard linear stability analysis of the 
above equations, which is analogous to the analysis for 
a monophase front with the added ingredient that the 
Gibbs-Thomson condition is coupled to a slow evolution 
equation for X(x, t) obtained by combining Eqs. (90) and 

©■ 

We start the stability analysis by writing the pertur- 
bations £ and y in terms of Fourier modes, 



£(x, t) = £fc exp(ikx + tot) 
y(x, t) = yt exp(z/cx + tut). 



(92a) 
(92b) 



The impurity diffusion field is expanded in Fourier modes 
according to Eq. (p4), and carrying out the same steps as 
from Eqs. (p4j) to (|58|) allows us to determine the Fourier 
coefficients. a result, Eq. (Bq) is replaced by 



bk = 



qk - ~(1- k E ) 





4M 




I 2 J 



with 



1 



(93) 



(94) 



We will in the following again use the quasistationary 
approximation of the impurity diffusion equation, which 
corresponds to dropping the term oj/D on the RHS of 
Eq. (|94|). As discussed before, we are mainly interested 
in perturbation wavelengths much smaller than the dif- 
fusion length. Note that in this limit, and within the 
quasistationary approximation, we have qk ~ \k\. We 
will not, however, make use of this simplification for the 
sake of generality. Next, we linearize the JH formula 
around the initial spacing Aq and the pulling speed v p : 



AT JH (X,v) = AT m (X ,v p ) 
, dATm 



dX 
<9ATjh 



\ ,v, 



\ dy 
ox 



dv 



(95) 



\ ,v p 



where we have used v„ — v p + £ and A — Ao = X^dy/dx 
from Eq. ( |90| ) . The Gibbs-Thomson condition, Eq. 
linearized in the perturbations, becomes 
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FIG. 12. Comparison between the full calculation and the 
long- wavelength limit, Eq. (^), for the stability spectrum. 
Symbols: full calculation, Eq. (fr4|), circles: real modes, 
crosses: complex modes. Lines: long-wavelength limit, Eq. 
([)?]), solid lines: real modes, dashed lines: complex modes. 
The parameters are the same as for Figs. fj(f) and 



G£ k = mAc [ -j£k — bk j 7^— 



v p k 2 
« 



dAT m 



dv 



(96) 



Inserting Eq. (93) for bk with the quasistationary approx- 
imation for qk, and using Eqs. ( |27| ) and ( |28| ) to calculate 
the derivatives of AXjh, we obtain a quadratic equation 
for u>. In the dimensionless quantities defined as before, 
this equation reads 



Po(k) 2t7(1 — 77) 



fi" 2 



g 2wPer 2 k E 

wr ~ o — 1 n 

2 po(«) 



2n 2 P(r 1 ) 
? ? (1 - 77) 



1 

A2 



,.2 



0. 



(97) 



with po(«0 defined by Eq. (Q). Figure p"2| shows a com- 
parison between this formula and the full calculation, 
Eq. fl74|), for the same parameters as in Fig. ^[ We see 
that, indeed, all features of the spectrum at small k are 
correctly predicted, including the transition from real to 
complex growth rates with growing k. This means that 
the simple calculation outlined above is able to capture 
the qualitatively new feature of the instability. Formally, 
the occurrence of the complex growth rates is due to the 
fact that the equation is quadratic in u>, whereas the anal- 
ogous equation of the Mullins-Sekerka calculation is lin- 
ear. This difference arises from the growth constraint 
resulting from Cahn's hypothesis. Physically, the change 
in the local lamellar spacing resulting from this constraint 
counteracts the destabilization of the front by the impuri- 
ties. This is due to the fact that in a convex (i.e. protrud- 
ing) part of the front, the growing lamellar spacing leads 
to an increase in the JH undercooling, whereas the inverse 
is true for concave parts. The magnitude of this effect 



is proportional to the slope of the JH undercooling ver- 
sus spacing curve, which increases with lamellar spacing. 
The oscillations occur because only the time derivative 
of the spacing depends on the instantaneous front shape, 
but not the spacing itself. Therefore, the perturbation of 
the lamellar spacing is phase shifted by it/ 2 with respect 
to the perturbation of the front shape. Consequently, 
during an oscillation cycle the former has its maximum 
amplitude when the front is planar. In summary, the 
long-wavelength oscillations are created by the interplay 
of the destabilizing effect of the impurity diffusion field, 
which is the same as for a monophase solidification front, 
and the dynamical response of the underlying lamellar 
structure. 



B. Inclusion of surface tension 

As shown in Fig. ^|, the spectrum derived from the 
stability analysis of the effective interface free boundary 
problem is in good quantitative agreement with the full 
spectrum at small k. This approach, however, fails to 
predict the restabilization of the interface at larger k be- 
cause it lacks capillarity. To add this effect, we can use 
the insights of Sec. V for the symmetric case, where it 
was noted that the eutectic stability spectrum could be 
interpreted as a planar interface spectrum with an effec- 
tive surface tension. This suggests that we can simply 
add to the Gibbs-Thomson condition (Eq. |8S[ ) a capil- 
lary term proportional to the curvature of the effective 
interface, which yields the new condition 



T = T E - mc~ ATj H (A,u n ) - T cS K[£], 



(98) 



where T e g is an effective Gibbs-Thomson coefficient. For 
the completely symmetric case, we can identify F e g by 
comparing Eq. (|97|) to Eq. (75) with the approximate 
expressions Eqs. ( |76[ ) for the coefficients a(n), 6(k), and 
c(/c), which yields at once the expression 



MAcPcAp 

I cff — I E H ^5 



(99) 



where = Tcos^ and Rq is given by Eq. (Bll 



Let us now consider the extension of this result to a 
general alloy phase diagram and an arbitrary composi- 
tion. The RHS of Eq. (|9^) contains two contributions 
that arise from the large-scale bending of the effective 
interface. The first is, in terms of the discrete formal- 
ism, the part of the curvature matrix K that is not con- 
tained in the JH formula. For arbitrary Gibbs-Thomson 
constants, contact angles, and volume fractions, simple 
arguments detailed in appendix C lead to the conclusion 
that the reaction of a composite interface to a small cur- 
vature can be descr ibed by a Gibbs-Thomson constant 
Te defined by Eq. (C4) that depends on T a , Tp, 8 a , 8p 
and the volume fraction 77. 

The second term on the RHS of Eq. ( p9| ) originates 
from the eutectic diffusion field, which was shown in Sec. 
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FI G. 13. Comparison between the full calculation and Eq. 
(L04) for an off-eutectic composition (77 = 0.65). The other 
parameters are as for Figs, g and [l^. Symbols: full calcula- 



tion, Eq. (p 
Lines: Eq. 
plex modes. 



(circles: real modes, crosses: complex modes), 
solid lines: real modes, dashed lines: com- 



V to have a stabilizing effect analogous to a supplemen- 
tary capillary term. For general alloy composition and 
phase diagram, this contribution can in principle be ex- 
tracted by expanding the complete spectrum to order k 2 . 
Since, as noted earlier, calculating this spectrum involves 
finding the roots of a fourth order polynomial in fi, this 
expansion is extremely tedious and was not carried out 
here. We have found numerically, however, that reason- 
ably accurate predictions can be obtained for off-eutectic 
compositions if we sim ply use T c g defined by Eq. ( J9S| ) 
with given by Eq. (C4). In view of the large uncer- 
tainty in the knowledge of several of the materials param- 
eters, notably the Gibbs-Thomson constants and the dif- 
fusion coefficients, this level of accuracy seems presently 
sufficient to interpret experimental results. Only very 
precise experiments could probe the differences between 
the full calculation and this approximation. 

Let us state the final result for the stability spectrum in 
two different forms to display the analogies with the MS 
and DL calculations, respectively. First, to compare to 
the MS instability, the dispersion relation can be written 



In the dimensionless quantities used by DL, the result 



is 



wr 



P(ji) 



po(n) 277(1-77) 



9/ 



g f ^K 2r )(r\)P(r\) \ i<- 



wr , 

2 V 

2wPer 2 k E ~ 



A 2 



Bo — 



Po(«) 
2tt 2 P(t]) 
77(1 - 77) 



1 

A 2 



1 - io ) ^ = 0, 



with the function 7(77) given by 



7(?7) 



T E {ri){m a + mp)/2 



(1 — rj)mpT a sm6 a + 77TO Q r^sin( 



(104) 



(105) 



Figure |l3| sh ows both the full calculation and the result 
of Eq. (104) for the stability spectrum of the symmetric 
phase diagram at an off-eutectic composition, 7/ = 0.65. 
The two are in reasonably good quantitative agreement, 
even though the value of R Q corresponding to 77 = 1/2 
was used in Eq. (|99|). 

Let us now briefly indicate which modifications will oc- 
cur if the two impurity liquidus slopes differ, rh a ^ rhp. 
In this case, the eutectic composition depends on the 
impurity concentration, and there is a eutectic bound- 
ary layer with a magnitude depending on the impurity 
concentration at the interface. If the two diffusion co- 
efficients D and D are equal (r = 1), all results carry 
over if we use the liquidus slope M defined by Eq. (|5|) 
instead of fh in all equations, and redefine the parameter 
w = MAc/MAc. This should usually be a reasonable 
approximation. If the two diffusion lengths are very dif- 
ferent, however, one would have to consider the eutectic 
boundary layer separately. Then, a separate Fourier ex- 
pansion has to be used for the eutectic boundary layer, 
and Eq. (|9^) is replaced by a more complicated form 
containing both diffusion lengths. 



^l = A(kl)(l- — -^(kl) 2 
2D K '[ 2l T 2i y ' 



2fe, 



lA(kl) 
4?7iAc 



AT ■ 

t—i-L mir 



Ay uj 



u)A 

v p 



(100) 



where A(kl), the thermal length It and the effective cap- 
illary length do are defined by 



A(kl) = yi + {kl) 2 
l T = mAc/G, 

do = -r-rz- 

mAc 



l + 2fe 



E j 



(101) 
(102) 

(103) 



Without the eutectic part on the RHS, Eq. (100) is the 
classical MS result for the one-sided model. 



VII. SHORT- WAVELENGTH MODES 

Up to now, we have only considered the instabilities aris- 
ing from the diffusive and MS modes. Let us now turn 
to the effect of the ternary impurity on the 2A-oscillatory 
(2A-Q) instability at off-eutectic compositions. In Fig. 
0, we show a series of spectra at an off-eutectic compo- 
sition (77 = 0.8) with increasing impurity concentration. 
The first spectrum, without impurity, again reproduces 
one of DL's figures. The diffusive branch is completely 
stable, but there is an unstable complex branch, with the 
most unstable mode at k — 0.5. For a small impurity con- 
centration (w = 0.01), the long- wavelength morpholog- 
ical instability is simultaneously present, but the struc- 
ture of the spectrum stays qualitatively unchanged. For 
still higher impurity concentration (w = 0.05), however, 
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the spectrum becomes quite different: a singie branch 
of the spectrum now contains both the long- and short- 
wavelength instabilities, making their distinction some- 
what arbitrary, whereas the other branch is completely 
stable. We have always found a similar structure for im- 
purity concentrations larger than w « 0.02. Figure |l4|(d) 
shows only the unstable branch for g = 0.05, w = 0.1, 
and A = 1, for a comparison with the spectrum at the 
eutectic composition shown in Fig. ^|(b). The long- 
wavelength part of the two spectra is very similar, but at 
the off-eutectic composition the most unstable mode is 
the 2A-0 mode. The growth rates of the long- and short- 
wavelength instabilities, however, are not very different, 
and we can expect a competition between the two. 

At K = 0.5, the matrix elements of A become real. 
The characteristic equation can again be factored in two 
quadratic equations, which are simply 



A a ' a = and A 13 ' 13 = 0. 



(106) 



This allows us, in particular, to obtain an equation for 
the neutral stability boundaries where the 2A-0 mode 
first becomes unstable. For the model alloy with the 
symmetric phase diagram, we have on the a-rich side of 
the phase diagram (77 > 0.5) 



g = 2wr 2k E rPeS 1 (0.5, 77) + 2r? - 1 

77 



+ l{R 2 (r l )-S 3 (0.5 ) r l )-^-P(ri) 



(107) 
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where i?2(^) and S^k, 77) are defined in appendix A. This 
is a direct generalization of DL's result for the binary eu- 
tectic. Figure [l^ shows the resulting stability diagram 
for g — 0.2 and A = 1. The dashed line is the consti- 
tutional supercooling criterion, and the long-range insta- 
bility is present everywhere a bove this line. The solid 
line was calculated using Eq. ( |l07| ), and the 2A-0 mode 
is unstable to the right of this line. We see that when 
the impurity concentration increases, the range in vol- 
ume fraction for which the eutectic front is stable de- 
creases. This means that, not surprisingly, the ternary 
impurity boundary layer enhances the oscillatory insta- 
bility. Furthermore, there is a large region in parameter 
space where the two instabilities compete and the fastest 
growing linear mode needs to be identified from a plot of 
the spectrum. 



VIII. CONCLUSION 



FIG. 14. Stability spectra for 77 = 0.75, A = 1, 6 = 45°, 
k E = 0.05, and (a) g = 0, w = 0, (b) g = 0, w = 0.01, (c) 
g = 0, w = 0.05, and (d) g = 0.05, w = 0.1. In (d), only the 
unstable branch is shown. 



We have performed a linear stability analysis of a thin 
lamellar eutectic interface in the presence of a ternary 
impurity to investigate the initial stages of colony forma- 
tion. The extension of Datye and Langer's method has 
allowed us to calculate the complete stability spectrum 
of the steady-state interface. From previous numerical 
studies of the binary eutectic case 21 1, we expect this 
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FIG. 15. Stability diagram for the symmetric eutectic al- 
loy in the space (t),w) for A = 1, g = 0.2, Pe = 0.01, and 
fee = 0.05. Solid line: neutral stability limit for the 2A-0 
mode. Dashed line: constitutional supercooling criterion. 

discrete stability analysis to be quantitatively accurate 
for spacings near the JH minimum undercooling spacing 

Amin ■ 

The most dramatic conclusion resulting from our anal- 
ysis is that the morphological instability of the eutectic 
interface induced by the ternary impurity is oscillatory, 
in contrast to the standard MS instability of a planar in- 
terface for a dilute binary alloy, which is non-oscillatory. 
We have seen that oscillatory modes originate from the 
interplay between the diffusive instability driven by the 
ternary impurity and the 'dynamical feedback' of the lo- 
cal change in lamellar spacing on the front motion. In a 
transient regime, these modes should create oscillatory 
microstructures with a wavelength of several lamellar 
spacings such as the ones displayed in Fig. |l^. There in- 
deed seems to be recent experimental evidence for large- 
scale oscillatory structures of this type in a transparent 
organic with a dilute ternary impurity |33f| , but a more 
detailed comparison between theory and experiment is 
now needed. 

Aside from its oscillatory character, the morpholog- 
ical instability of the eutectic interface is qualitatively 
similar, near onset, to the standard MS instability of a 
monophase front. In particular, we find that the expres- 
sions for the critical onset velocity and morphological in- 
stability wavelength are analogous to those for the clas- 
sic Mullins-Sekerka instability of a planar interface. The 
main difference is that the restabilization of the interface 
at short wavelength is controlled by an effective surface 
tension that depends on the geometry of the lamellar in- 
terface and on intcrlamcllar diffusion, which has a resta- 
bilizing effect. One consequence of this result is that the 
constitutional supercooling criterion that has been com- 
monly used in the metallurgical literature to predict the 
onset of instability is indeed applicable for typical alloy 
compositions. Note, however, that this criterion becomes 
inaccurate for very small concentrations of the ternary 
impurity. 



Above the onset of instability, the stability spectrum 
can exhibit both real and complex modes. The scaling of 
the wavelength of the fastest growing mode with pulling 
velocity depends on the nature of the mode. For com- 
plex modes, far above the onset this wavelength scales 
as the geometric mean of the capillary length and the 
diffusion length. For real modes, the situation is more 
complicated. In both cases, for fixed velocity and lamel- 
lar spacing, the wavelength far above the onset and at 
sufficiently high impurity concentration scales as the in- 
verse square root of the impurity concentration. Note 
that all these statements concern the primary instability 
of the eutectic front and not the finally selected colony 
spacing. After the first colonies have formed, the non- 
planar front may undergo a complicated sequence of cell 
elimination or tip splitting events, as during the develop- 
ment of monophase cellular structures. 

Furthermore, we have shown that the eutectic front 
dynamics on scales much larger than the lamellar spac- 
ing can be formulated as a free boundary problem with 
a modified Gibbs-Thomson condition that is coupled to 
a slow evolution equation for the lamellar spacing. This 
formulation provides a deeper physical understanding of 
the eutectic front dynamics on this scale. In addition, 
we have shown that it can be used to calculate an ap- 
proximate stability spectrum that is well suited to in- 
terpret experimental data. The effective capillary length 
appearing in this spectrum contains contributions both 
from an averaging over the material properties of the 
two phases, weighted by geometric factors, and from the 
eutectic interlamellar diffusion field, which acts as a sta- 
bilizing force. 

Finally, we have found that the short-wavelength os- 
cillatory instability, already present in a binary eutec- 
tic, is enhanced by the ternary impurity boundary layer. 
This reduces the composition range for stable lamellar 
eutectic growth, even below constitutional supercooling. 
Above constitutional supercooling and for sufficiently off- 
eutectic compositions, the long-range and 2A-oscillatory 
instabilities are both present and may compete with each 
other. 

In conclusion, we have shown that the instability of a 
lamellar eutectic interface in the presence of a ternary im- 
purity is in some respects similar to the Mullins-Sekerka 
instability of a dilute binary alloy, but also presents 
striking differences. There are two interesting future 
prospects. First, dynamical simulations of the complete 
equations of motion are necessary to go beyond this linear 
stability analysis and to investigate the subsequent stages 
of the instability, as well as to determine what structures 
are ultimately formed. Work on these issues using the 
phase-field method is currently in progress. Second, it 
seems worthwhile to extend the effective interface ap- 
proach to a nonlinear regime to model the shape and 
dynamics of fully developed colonies, as depicted in Fig. 

B 
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APPENDIX A: SUMMARY OF DATYE AND 
LANGER'S RESULTS 

We will state here DL's results for the matrices G, K, and 
U and transform them into our notations. The matrix 
G can be simply read off the definitions of the average 
interface positions (Qj, Eqs. (p8|), and of the Fourier 
expansion Eq. (|4l|): 



G = 



G 



1 1 

,ik\ 1 



(Al) 



The curvature matrix contains two contributions. The 
first arises from the change in the local lamellar spac- 
ing due to the horizontal displacements The second 
appears when the interface is bent on a scale of sev- 
eral lamellae. Then, the trij unction points are turned 
by small angles with respect to their steady-state orien- 
tation. These angles can be related to the time deriva- 
tives of the horizontal displacements. Finally, we use the 
growth constraints Eq. (MG) and obtain 



2ie lkx ' 2 sin(fcA/2) 

a 2 

_ 2ie ifcA / 2 sin(fcA/2) 
" A 2 



2v p sin ( 



cos 6 n 



(A2) 



2v p sin 6p 
wA(l - rj) 2 



COS ( 



(A3) 



Note that in the general case the scaling with respect to 
the physical parameters would remain the same; however, 
additional coefficients depending on the angles 6 a and 9p 
and the liquidus slopes m a and mp would appear. 

For the matrix U, we will just state DL's results; for 
more details, see j22j. The calculation is straightforward 
but tedious because we have to treat the interlamellar 
diffusion. This brings in various sums over the Fourier 
modes of the steady-state expansion, Eq. (121). We define 
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where the asterisks in the last two equations denote com- 
plex conjugation of all the coefficients, but not u). 

For the symmetric phase diagram, we may conve- 
niently rewrite this matrix in the dimensionless parame- 
ters defined in Sec. IV. We remark that in this case we 
have 



IT sin 9 
MAcP(rj)' 



(A6) 



This relation can be used to eliminate sin# in the ma- 
trix; the reduced lamellar spacing A appears, and, for 
example, we obtain for K°' a 



K 
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(A7) 
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Using these quantities, the matrix elements of U are 
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VM = -±(ilU?{K,ri) + Ui{K,ri) 



+ — — sin{TTK)U 3 (K,l-v) ), (A18) 



= \J a - a * (A19) 



with 



r?) = - [A a Si(n, rj) - A p S^K,r,)] (A21) 



U l( K iV) = i A P S l( K ' i-V)- A a S%{K, 1 - 7?)] 

(A22) 



tf 2 a («,»?) = -S + - [SZ(k, V ) - R 2 (ri)] (A23) 



Ug{n, v) = B Q + ^— y 1 - r?) - R 2 (l - t,)} 



(A24) 



U 3 {K,r]) = - [2P( V )/r ] -R 1 (r 1 ) - 2SZ(k, V )] . (A25) 



APPENDIX B: LIMIT OF THE DL SUMS FOR 

K — ► 

For the detailed study of the symmetric phase diagram at 
eutectic composition, we need to know the leading order 
behavior of the functions U°, U% , and Us, for k — > 0. 
To this end, we have to expand the sums S n in k and 
resum the resulting terms. We will systematically neglect 
terms of relative magnitude Pe, as was already done in 
DL's calculations leading to the results of the preceding 
appendix. 

We have to single out terms containing the function 
Po( K ) m the denominator, because these terms will be 
singular in the limit k — > 0. To see this, note that we 
have 

, , f 27r 2 K 2 /Pe + C(k 4 ) for 2ttk < Pc 
P°W = \2*\k\ for 2™»Pe ( B1 ) 

We will be interested in a regime where the wavelength of 
the perturbation is larger than the lamellar spacing, but 
much smaller than the diffusion length; hence the latter 
limit applies. 



Similarly, the function po{n) will become small when 
k tends to 0: 

-,>._/ rPek E + 27r 2 K 2 /Pe + 0{k a ) 2ttk < Pc 

(B2) 

Hence the terms proportional to p^ 1 in the impurity con- 
tributions have to be considered separately. 
Expanding S\{K,rj) and ^(fi;, 77), we obtain 
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The expansions for the impurity sums Si(K,rj) and 
52(^,77) are obtained by replacing po(n) by Po(k) in the 
above expressions. 

With the help of these expressions and the definition 
of [7f, we find for 77 = 1/2 and A a = Ap = 1/2: 



Rc(e- 47rK / 2 J7f( K ,l/2)) 
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Rather remarkably, all the singular terms cancel out. A 
comparison with direct numerical summation shows that 
the limit behavior is correct and that the neglected terms 
sum up to a correction that does not exceed P(l/2), even 
for large values of k. 

For the impurities, we need the expression 

§! (k, 1/2) cos ^ + Re (e- l7TK ' 2 S* (k, 1/2) 



1 - ^7T 2 K 2 



1 



9 - 7 x + ^ [^(V2) + #3(1/2)] ttV 
+ O(k 4 /p ) + O(k 4 ). 



(B9) 
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We will simplify our task for the function U%, which 
contains the most difficult sum S3 (re, 77), by directly ex- 
panding the product e I7rK / 2 S3 (re, 1/2). Using P = for 
the symmetric phase diagram at eutectic composition, we 
obtain 



Re 



(e^'Ujj = P re 2 /2 + G(re 4 ) 



(BIO) 



with 



#0 = 2tH 



3^(1/2) - -R 2 {l/2) - P 7 (l/2) 



0.4965 



2 sin(7r?7m) 
(7rm) 2 



(7m) 4 



(Bll) 



(B12) 



— sin(7T7yn) 
n 



cos 2 (7r(n — m)/2) 
7r 2 (n — to) 2 



(B13) 



Comparison to direct summation shows that the expres- 
sion ( B10| ) is accurate to within 5% over the whole range 
of re. 

Finally, to express U3, we need the expansion 



S-4(re, 77) 



rj + inr]' 2 K 2-K l j] i K 2 + in 3 r] 4 K, 3 



Po(re) 3p (re) 
+ Ri(v)/2 + 2m [P(n) - vRiiv)} « 
+ [6i? 4 (??) - 877^3(77) - 77^(77)] 7r 2 re 2 /4 
+ C(re 3 ). (B14) 

ForEq. @, we need 

2sin(7rre)Re (ie l7rK/2 [/ 3 (re, 1/2)) = 87t 2 P(?7)k 2 + O(re) 4 . 

(B15) 

Collecting all these results, we can finally write down 
the complete expressions for the coefficients a(re), 6(re), 
and c(re) of Eq. (|75|) up to order re 2 : 



a(re) =2P(l/2) + 
1 - 



12J2b(1/2) - ^(1/2) 



2 2 
7r re 



24^ K 



po( K ) 



6(re) 



+ [P(l/2) + P 3 (l/2)]^ 2 re 2 

(tor- g/2) (1 -7t 2 k 2 /8) 
'27r 2 P(l/2)cot6» P 
A 2 + T 



(B16) 
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c(re) = 8tt^P(1/2) 1 



Po(re) 

[P(l/2) + P 3 (l/2)K 2 re 2 
1 

A 2 



(B17) 
(B18) 



To obtain the expressions given by Eqs. (|7q), we remark 
that for small re we can neglect most of the terms listed 
above. We have to be careful, however, to keep track of 
all physical effects. For example, formally the leading or- 
der term in a(re) is of order re -1 if we use po(re) = 27r|re|. 
But this term comes with a prefactor w, proportional 
to the impurity concentration, whereas the leading order 
term in the expansion of U", arising from the eutectic 
diffusion field, is independent of w. Hence for the ap- 
proximation to be valid for a rbitr ary w, we need to keep 
both terms, leading to Eq. (76a) for a(re). Similarly, in 
6(re), we need to keep the dominant terms for each phys- 
ical effect, even if their order is higher than other terms 
we may neglect. In particular, we must keep the capil- 
lary term that is of order re 2 and has a prefactor of order 
unity. Let us show that, at the onset of instability, we 
can indeed neglect other terms of order re and re 2 . Keep in 
mind that we are interested in a regime where re is small, 
but not too small; a plausible estimate is re « 0.01. First, 
there is the correction —n 2 K 2 (wr—g/2)/8 to the constitu- 
tional supercooling criterion, of order re 2 . But as wr—g/2 
is very small near the onset, this term is actually much 
smaller than the capillary term. Next, there are terms 
of orders re and re 2 arising from the impurity contribu- 
tion. But all of these come with a prefactor of wPekE', as 
both w and Pe are small quantities and fee < 1 , it seems 
justified to neglect them. For example, for w = 0.1, 
Pe = 0.01, and re = 0.01, the largest neglected term is 
of order wPere ~ 10 -5 , wherea s the capillary term is of 
order re 2 ~ 10~ 4 . Hence, Eq. ( [76b] ) for 6(re) seems well 
justified at the onset of instability. 



APPENDIX C: EFFECTIVE SURFACE TENSION 

We will give here an expression for the Gibbs-Thomson 
constant of a lamellar eutectic interface T e that describes 
the shift of the average interface temperature when the 
composite interface is curved on a scale much larger than 
the lamellar spacing. This analysis is necessary because, 
in a composite material, the interface with the weaker 
surface tension will absorb more of the curvature, leading 
to an effective surface tension that depends on the volume 
fraction. Note that the expression derived here is valid in 
thermodynamic equilibrium and contains only the "ge- 
ometric part" of the effective Gibbs-Thomson constant 
r e ff for a moving eutectic front, in which the stabilizing 
effect of the intcrlamcllar diffusion has to be included. 

Consider a lamellar interface that is curved such that 
the a(3 solid-solid interfaces on the two sides of a lamella 
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pair make a small angle <f>. Suppose that the [3a interface 
between them is turned by an angle <pi . The corrections 
of the curvature with respect to the planar front values 
are then given by (using the fact that <j) ~ A/i2 for large 
radii of curvature R) 



6K n 



cos 9 a 



)i and 5Kp 



cos Op 



(<f>-<Pi). (ci) 



As the average temperature of neighboring lamellae 
should be the same, we must have 



5K a T a = SKpTp. 



(C2) 



From this condition, we can determine the unknown an- 
gle 4>\. Finally, we obtain the undercooling of the inter- 
face as 



with 



AT = T E K 



r a T(3 COS 9 a COS 9/3 

(1 - r])T a cos9 a + rjTp 



COS 0/3 



(C3) 



(C4) 
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